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Abstract 



The straightforward supersymmetrization of the Standard Model (SM) results in a phe- 
nomenologically inconsistent theory in which Baryon number (B) and Lepton number (L) 
are violated by dimension 4 operators, inducing fast proton decay. Proton stability allows 
only for separate L or B violation and, if neutrinos arc massive Majorana particles, L 
violating terms must be present. In this thesis I will study a Supersymmetric Standard 
Model (SSM) realization with B conservation and minimal L violation. In this frame- 
work L is mildly violated only by super-renormalizable terms, allowing for small neutrino 
Majorana masses. This model is more predictive than the Baryon-Parity SSM. The in- 
duced dimension 4 L violating couplings are not arbitrary, and automatically satisfy all 
experimental constraints. After introducing the theoretical framework for supersymmet- 
ric models without Lepton number, I will discuss the phenomenology of the (unstable) 
lightest neutrahno and of the hghtest stop. I will show that the leptonic decays of the 
stop can be related to the neutrino parameters, and in particular their measurement can 
indirectly probe the size of the solar neutrinos mixing angle. 
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Resumen 



El Modelo Estandar Supersimetrico (MES) es el modelo supersimetrico con mmimo 
niimero de particulas correspondientes a aquellas del Modelo Estandar (ME), un doblete 
de Higgs adicional y todos los companeros supersimetricos de las particulas del ME. El 
MES permite la aparicion simultanea de tcrminos que violan numero leptonico y barionico 
y de este mode es excluido debido a que el proton decae a traves de interacciones elec- 
trodebiles. Sin embargo es bien conocido que la estabilidad del proton permite la presencia 
ya sea de violacion de numero leptonico o de la violacion de numero barionico. De he- 
cho, si los neutrinos tienen masa como es sugerido por las anomali'as solar y atmosferica 
no existe una razon de peso para asumir que los terminos que violan numero leptonico 
esten ausentes. En este trabajo estudiaremos la realizacion del MES con conservacion 
de numero barionico y violacion minima de numero leptonico. Los linicos terminos con 
violacion de numero leptonico en este caso provienen del termino de masa mas general 
que permite la invarianza gauge en el superpotencial del MES. Al modelo correspondiente 
le llamaremos Modelo Estandar Supersimetrico Superrenormalizable (MESS). El MESS 
es mas predictivo y teoricamente mas atractivo que el MES con conservacion de numero 
barionico pero con terminos arbitrarios de violacion de numero leptonico y en general 
satisface todos las restricciones que surgen de esta violacion. 

En el MESS la Particula Supersimetrica mas Liviana (PSL) puede decaer. En par- 
ticular estudiaremos la fenomenologi'a esperada del neutralino en LEP2. En el caso de 
particulas supersimetricas con mayor masa que la PSL, los decaimientos que violan numero 
leptonico pueden Uegar a ser detectables en los aceleradores, incluso en el caso de masas 
de neutrinos tan pequenas como las sugeridas por las anomalias de neutrinos. Especial- 
mente interesante en este caso son los decaimientos que violan numero leptonico del stop 
(companero supersimetrico del quark top), y que seran estudiados en este trabajo. Estos 
decaimientos implican la posibilidad de probar el angulo de mezcla solar y asi la posibil- 
idad de relacionar la fisica de colisionadores de alta energia con los experimentos sobre 
oscilacion de neutrinos solares. 
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Introduction 



The Supersymmetric Standard Model (SSM) [1, 2] is the Supersymmetric (SUSY) model 
with minimal possible number of particles, corresponding to those of the Standard Model 
(SM), one extra Higgs doublet and all the superpartners of the SM particles. We define 
the Minimal Supersymmetric Standard Model (MSSM) [2, 3] as the SSM with minimum 
possible number of couplings. 

Unlike the MSSM, the Supersymmetric Standard Model allows many gauge invariant 
terms violating baryon number (B) and lepton number (L). Consequently, this most 
general case is excluded because the proton would decay with a weak decay rate. However, 
it is well known that the stability of proton allows either L ot B violation, because the 
nucleon decays normally requires not only the baryon number violation but also lepton 
number violation. In the usual way to avoid proton decay, all such terms are forbidden by 
imposing some X symmetry [4]. For example the anomaly free gauge discrete symmetry 
known as /^-parity [5-7], forbids all the renormalizable L and B violating terms. It is 
worth to stress, in that case it is possible have L violation at the low energy theory. 
For example, superrenormalizable (bilinear) L violating terms in the low energy SSM 
may appear from renormalizable (trilinear) terms allowed at some high energy where the 
i?-parity is still unbroken [8]. 

From a phcnomcnological point of view, it is most important to ensure that there 
are no interaction terms in the lagrangian which lead to rapid proton decay and, in this 
respect, other discrete symmetries can be used which are even more effective than im- 
parity. For example the anomaly free discrete symmetries equivalent to baryon parity 
[7, 9-12] or lepton parity [10, 13]. The former forbids dimension-4 and 5 B violating 
operators while the latter forbids dimension-2, 4 and 5 L violating operators. In this way, 
general models of L or S violation are expected to come from X symmetries different 
from R parity. 

If neutrinos are massive Majorana particles, lepton number is violated, and there is 
no compelling reason to assume that L violating terms are absent from the superpoten- 
tial. Indeed, there is considerable theoretical and phenomenological interest in studying 
possible implications of alternative SUSY scenarios in which L is broken [14-25]. This is 
especially so considering the fact that it provides an appealing joint explanation of the 
solar and atmospheric neutrino [26] anomalies which has, in addition, the virtue of being 
testable at present and future accelerators like LEP [27-29], Tevatron[27, 30-33], LHC [25, 
34] or a linear e+e" collider [35]. The effects of L violation can be large enough to be 
experimentally observable. 

An special case of baryon parity SSM only contain dimension-2 L violating terms in the 
superpotential. They have mass coefficients protected by the non-renormalizable theorem 
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in SUSY [36] and are called superrenormalizable terms. In fact, once mass parameters 
are introduced, the term 

Wf, = HaLaHu (1) 

(rather than only fioHdHu) is the most natural choice. Where Lq = and a = 0, 1, 2, 3. 

The L violation through only superrenormalizable terms could arise explicitly as in [8, 
16, 37-50] as a residual effect of some larger theory. Most of these models do not introduce 
new particles to the SSM and we will call them Superrenormalizable Super symmetric 
Standard Models (SSSM). 

Alternatively the L violation could arise spontaneously, through nonzero vacuum ex- 
pectation values (vev's) of singlet fields [51-57] which break i?-parity at low energies. 
Consequently these models contain additional fields not present in the SSM. Although 
in this work we will concentrate in the SSSM, most phenomenological features of the 
spontaneous L violating SUSY models are reproduced also for the SSSM 

SSSM is more predictive and theoretically more attractive than L violating SSM 
through arbitrary dimension-4 operators. Moreover, unlike this baryon parity SSM, the 
induced dimension-4 L violating couplings in the SSSM are not arbitrary and in general 
automatically satisfy all experimental constraints on L violation. This renders a system- 
atic way to study L violating signals [37, 58-64] and leads to effects that can be large 
enough to be experimentally observable, even in the case where neutrino masses are as 
small as indicated by the simplest interpretation of solar and atmospheric neutrino data 
[34, and references therein]. Moreover, the SSSM follow a specific pattern which can be 
easily characterized. These features have been exploited in order to describe the expected 
SSSM signals (see section 2.4) 

In chapter 1 we will study the SSM model for one generation. This "toy model" turn 
to be very illustrative in the discussion of key ingredients of the general baryon parity 
SSM, such as the definition of basis independent parameters and the generation of the 
tree level neutrino mass. 

In chapter 2 we emphasize the problems with the SSM and study in a systematic 
way the several alternative SSM models. Instead of the usual approach of ad-hoc matter 
parities we use theoretical best motivated anomaly free gauge discrete symmetries. They 
allow an easier connection with other realistic ways of generate the unknown X symmetry 
responsible for the proton stability in the SSM, such us R symmetries, gauge symmetries, 
flavor symmetries, Peccei-Quinn symmetries. Next we review how the SSSM can be 
generated with all these X symmetry possibihties. The SSSM is then presented as the 
best motived and minimal realization of the SSM when the experimental evidence on 
neutrino masses is taken into account. 

In chapter 3 we build the mass spectrum of the SUSY particles in the SSSM. In the 
next chapters we study the phenomenology of the neutralino and the stops in the SSSM. 
The neutralino and the two body stop decay studies were performed in the supergravity 
version of the SSSM, while the three body decay of the stop was studied with arbitrary 
low energy parameters as inputs. 

In the neutralino case we will present a detailed study of the Lightest Supersymmetric 
Particle (LSP) decay properties and general features of the corresponding signals expected 
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at LEP2. It is well known that in models with Gauge Mediated Supersymmetry Breaking 
(GMSB) the lightest neutralino decays [65, 66], because in this case the gravitino is the 
LSP. We therefore also discuss the possibilities to distinguish between GMSB and the 
SSSM. 

We will study the decay modes of the lightest top squark in supergravity models where 
supersymmetry is realized with the SSSM. In such models the lightest stop could even be 
the lightest supersymmetric particle and be produced at LEP or Tevatron. Neither e^e~ 
collider data [27-29, 67] nor pp data from the Tevatron [68, 69] preclude this possibility. 
In contrast with ref. [70] here we focus in the effective model where the L violation is 
introduced through an explicit superrenormalizable term. This is substantially simpler 
than the full majoron version of the model considered previously. 

In order to discuss stop decays we also refine the work presented in Ref. [71-76] by 
giving, for the first time, an exact numerical calculation for the Flavor Changing Neutral 
Current (FCNC) process i cx?. We also compare the results obtained this way with 
those one gets by adopting the usual one-step or leading logarithm approximation in the 
Renormalization Group Equations (RGE). In contrast with the MSSM such an approxi- 
mation would be rather poor for our purposes, since we will be interested in comparing 
FCNC with L violating stop decay modes. Moreover, in contrast to ref. [70], where the 
magnitude of the stop - charm - neutralino coupling was a phenomenological parameter, 
here we assume a minimal supergravity scheme with universality of soft terms at the 
unification scale in which this coupling is induced radiatively and thus calculable. As we 
will see this has important phenomenological implications, for example in the behavior 
of the stop decays in the SSSM with respect to tan /3. We calculate its magnitude using 
a set of RGB's in which the running of the Yukawa couplings and soft breaking terms 
is taken into account. Here we also provide the analysis of the relationship of the stop 
decays in the SSSM with the magnitude of the heaviest neutrino mass. Motivated by the 
simplest oscillation interpretation of the Super-Kamiokande atmospheric neutrino data, 
we also generalize the treatment of the L violating decays by explicitly considering the 
case of light v^, masses, not previously discussed. 

In the SSSM, the stop can have new decay modes such as 

ti^hr (2) 

due to mixing between charged leptons and charginos. We show that this decay may be 
dominant or at least comparable to the ordinary L conserving mode 

ii^cxl (3) 

where denotes the lightest neutralino. 

Owing to the large top Yukawa coupling the stops have a quite different phenomenology 
compared to those of the first two generations of up-type squarks (see e.g. [77, 78] and 
references therein). The large Yukawa coupling implies a large mixing between it and Ir 
[79] and large couplings to the higgsino components of neutralinos and charginos. The 
large top quark mass also implies the existence of scenarios where all MSSM two-body 
decay modes of ii are kinematically forbidden at the tree-level (e.g. ti — > tx^bx^ ,tg). 
In such case higher order decays of ii become relevant [71, 73, 74, 80]: ii — > cxi^2-: 
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tl W^hx\^ tl H'^bxi, tl bl^ui, tl buil^, where / denotes e,yU, r. Also 
4-body decays may become important if the 3-body decays are kinematically forbidden 
[81]. In [73, 74, 80, 82] it has been shown that in the MSSM the three-body decay modes 
are in general much more important than the two body FCNC decay mode. Recently it 
has been demonstrated that not only LSP decays but also the light stop can be a good 
candidate for observing L violation, even if its magnitude is as small as indicated by the 
solutions to the present neutrino anomalies [31, 32, 70, 83]. In particular in [31] (see 
section 5.1.3) it has been demonstrated that there exists a large parameter region where 
the SSSM decay 

ii br 

is much more important than the MSSM decays 

It is therefore natural to ask if there exist scenarios where the decay ii — > 6 r is as 
important as the three-body decays. Note that in the SSSM the neutral (charged) Higgs- 
bosons mix with the neutral (charged) sleptons [84, 85]. These states are denoted by 
S^, Pj, and 5^ for the neutral scalars, pseudoscalars and charged scalars, respectively. 
Therefore in the SSSM one has the following three-body decay modes: 

ii ^ W+bx'i 

ii - Stbxl 
tl S^b ui 
ii bS°l+, 
tl ^ 6P°/+. 

We will show that there exist regions in parameter space where fi — > 6 t+ is sizeable and 
even the most important decay mode. In particular we will consider a mass range of 
tl, where it is difficult for the LHC to discover the light stop within the MSSM due to 
the large top background [86]. Contrary to the studies on LSP decay in the SSSM, the 
processes here are very sensitive to the value of the heaviest neutrino mass 



Chapter 1 



One generation Supersymmetric 
Standard Model 



In the Standard Model (SM) the radiative corrections to the Higgs mass are naturally 
of order of the Planck mass scale. However if there is a symmetry relating bosons and 
fermions, called supersymmetry, such as corrections turn to be at a controllable level. This 
symmetry at least doubles the field content of the SM. In this work we will study SUSY 
models with minimal possible number of particles, corresponding to those the standard 
model, one extra Higgs doublet and all the supcrpartners of the Standard Model (SM) 
particles. The most general renormalizable SUSY model respecting the gauge symmetries 
and with minimum field content will be called Supersymmetric Standard Model (SSM). 

The field content of the SM together with the requirement of Gsm = SU{3)c x 
SU{2)l X U{1)y gauge invariancc implies that the most general Lagrangian is character- 
ized by additional accidental U{1) symmetries implying Baryon number {B) and Lepton 
number (L) conservation at the renormalizable level. When the SM is supersymmetrized, 
this nice feature is lost. The introduction of the sTipcrpartncrs allows for several new 
Lorentz invariant couplings which are in confiict with present bounds in B and/or L 
violation. However in the one generation SSM B is automatically conserved (at the per- 
turbative level) while the L violation may be compatible with experimental constraints. 
In this chapter we study the one generation SSM and postpone the discussion of problems 
with the three-generation SSM for the next chapter. 



1.1 Supersymmetric Lagrangian 



A supersymmetric transformation in a realistic supersymmetric model, turns a bosonic 
state 4>i into its superpartner Weyl fermion ■0^ state and vice versa. It also converts the 
gauge boson field A" into a two component Weyl fermion gaugino A" and vice versa 
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1.1 Supersymmetric Lagrangian 



50i = eipi (1.1a) 

S{iJi)a = i{a^e\D^<Pi + e^Fi (1.1b) 

6F, = le^a^D^^, + V2g{T^(t^),e^ X^'' (1.1c) 

M» = -i=(eta^A» + At%6) (l.ld) 

= (^'^^"^e), F;, + 6„D« (Lie) 

SD" = ^ (e^a^^I^^A - D^X^&f'e) (l.lf) 



where the undotted (dotted in greek indices are used for the two components of 
the left-handed (right handed) Weyl spinors^. The are 2x2 matrices with Uo = (Ji 
being the identity and ai = —cJi the Pauh matrices. The index i runs over the gauge and 
flavor indices of the fermions (it is raised or lowered by hermitian conjugation); is an 
infinitesimal, anti commuting two-component Weyl fermion object which parameterizes 
the supersymmetric transformation. and are complex auxiliary fields which do not 
propagate and can be eliminated using their classical equations of motion. The index 
a runs over the adjoint representation of the gauge group under which all the chiral 
fields transform in a representation with hermitian matrices satisfying \T°', T^] = if°-^'^T'^. 
Finally the gauge transformations are 



D^(t>i^d^(t>i + igAl{T''(t>)^ (1.2a) 

D^i^i^d^i^i + igAliT'^i^), (1.2b) 

D^X^^d^X^-gP^'^AlX'^ (1.2c) 

f;, = d,Al - d^Al - gr'^A^Al (1.2d) 



where g is the gauge coupling. The gauge quantum numbers of ipi are the same as that 
its scalar superpartner (p^. For a complete discussion see [36]. 

As a result, in a renormalizable supersymmetric field theory, the interactions and 
masses of all particles are determined just by their gauge transformation properties and 
by the superpotential W (see for example [36]) 

W = \m'%'^^ + \y''%^j^k (1.3) 

where the superfield is a singlet field which contains as components all of the bosonic, 
fermionic and auxiliary field within the corresponding super multiplet, e.g. 0i D (0j, Fi) . 
W determines the scalar interactions of the theory as well the fermion masses and the 
Yukawa couplings. 

The superpotential, together with the SU{2i)c x SU{2)l x U{1)y gauge symmetry, 
lead to the following generic lagrangian 

^In general, a = 1, 2, . . . ,d = 2^^ whit N the number of supersymmetries. 
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-^SUSY — — 




where $i and are the SU{2)l doublets (see eq. (1.9)). 

Any reahstic phenomenological model must contain supersymmetry breaking. We 
use the usual approach of parameterized our ignorance of the specific mechanism by 
just introducing extra terms which break supersymmetry explicitly in the effective SUSY 
lagrangian. The additional possible soft supersymmetry terms in the previous lagrangian, 
assuming SU{3)c x SU{2)l x U{1)y gauge symmetry is 



1.2 One generation Supersymmetric Standard Model 

Since the third generation fermions have the heaviest masses in the Standard Model and 
the Lr violating processes are less restricted, it is often useful to make the approximation 
of keeping only the third family. 

The most general renormalizable superpotential for one generation of Standard Model 
quarks and leptons, imposing only gauge invariance is [41, 62, 63, 87-95] 



where Lq = Hs,. Note that there is no possibility of B violation in the one generation case. 
Moreover we have chosen a basis where the only allowed dimcnsion-4 L violating term, 
A333L3Q3D3, is absent from the superpotential. This minimizes the number of parame- 
ters in the superpotential because in converse case, the /is terms should reappear when we 
evolve the SUSY parameters to a different scale. The //q and //3 are mass terms that are 
protected by the supersymmetric non-renormalization theorem (for a review see [36]) . In 
particular, it means that once we have a theory which can explain why jia (« = 0, 1, 2, 3) 
is of order 10^ or lO^GeV at tree-level, we do not have to worry about fia made very 
large by radiative corrections involving the masses of some very heavy unknown particles. 



Koft = - 2 (^1^'-^' + + M3A3A3 + h.c) - m2.0*^0, 



(1.5) 




W = htQzU2.Hu + X0Q3D3L0 + hrLsEaLo - Aio^oi^2 - Ai3^3-f^2 



(1.6) 
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1.2 One generation Supersymmetric Standard Model 



Thus, the ^„ are called superrenormalizable terms to differentiate from the bilinear soft 
mass terms which are not protected by the supersymmetric non-renormalization theo- 
rem. We call the models where the L violation is induced only by Hi terms and have 
minimum content of fields {i — 1,2,3), Superrenormalizable Supersymmetric Standard 
Models (SSSM). In particular, in the one generation case, the SSM is equivalent to the 
SSSM. 

The soft potential in this case is 



^oft = £^ 



ah 



mass terms. (1.7) 



In this simple case it is easier the study of explicit calculations of L violating process 
(section 1.2.1), basis independent L violating parameters (section 1.2.2), and tree level 
neutrino mass generation (section 1.2.3). 



1.2.1 One generation SSM lagrangian 

The superpotential in eq. (1.6), can be rewritten as 

where 



Q3 - (tL bi) ^3 - (jTs fij 

which is of the form of the generic Lagrangian in eq. (1.4) if we define 



(01, ^1) 


= {h,tL) 


(05,^5) 




(09, ^9) 




(02,^2) 


= ih^b^) 


(06,V^6) 


= (-^«, H^) 


(0(10), V'(IO)) 


= (t'LiTl) 


(03,V'3) 




(07, ^7) 


= (-^0, -^0) 


(0(11), V'ai)) 


= i^W ^r) 


(04, -04) 


-ibR,bR) 


(08, -08) 


= (-^0 , -^0 ) 







and with this notation we have 

if i is odd. 



^0i+ly 





if i is even 




if i is odd, 
if i is even 



:i.8) 



:i.9) 



Comparing eqs. (1.3) and (1.8) we obtain the non-vanishing couplings 

In Appendix B.l we show an explicit calculation by using the previous formulas. 
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1.2.2 Basis independent petrameters 

In the low energy superpotential of the one generation SSM in eq. (1.6), it is possible 
to rotate away the L^H^ term from the superpotential by redefinition of and L3. It 
is worth stressing, however, that such a redefinition does not leave the full lagrangian 
(including soft breaking terms) invariant [16, 31, 37, 88, 89, 92, 95-101]. This generates a 
term and also affects the soft supersymmetry breaking terms. While the iisL^H^ 
term in the superpotential can be rotated away at a fixed energy scale, the corresponding 
L^Hu would still be present in the low energy theory. Alternatively, if one chooses to 
remove the fi^ term at the Planck scale, the terms L^H,^ will be radiatively generated due 
to the trilinear L violating term [16]. We illustrate this statement below in the context 
of the neutrino mass generation. 

The two bihnear L violating terms in eqs. (1.6) and (1.7) give rise to two different tree 
level contributions to neutrino mass: (1) The term /i^L^Hu provides a mixing between 
Ur and Hu- (2) The term BifiiL^Hu lead to a sncutrino vcv v?, which induce a mixing 
between z/^ and the gauginos. Note that the contribution to neutrino mass coming from 
the dimension-4 L violating term was rotated away. This is no longer possible in the full 
three-generation SSM and may lead to too large tree level neutrino mass [96] . 

The question is if all the bihnear terms can be rotated away from the SSM lagrangian 
so that, for example, the neutrino mass can arise only from the loop- level contribution 
from the dimension-4 L violating term. We will see that this can only be done under 
very specific conditions on the relevant soft parameters at the electroweak scale which 
should look like unacceptable fine-tuning. These conditions can be determined from the 
minimization equations with vs — in the basis where already the superrenormalizable 
term have been rotated away, inducing one trilinear violating term in the superpotential 
with coupling A3 = A333 . Moreover, the conditions turn to be not scale invariant because 
the Renormalization Group Equation (RGE) for /X3 have a contribution proportional to 
/^o [16, 88, 92, 96-99] 

167r2^ = 3/X0A0A3 (1.10) 
at 

so that, even when is zero to one scale, it is radiatively generated at another different 
scale. However if the is absent at some energy, this will never reappear in the 

superpotential because the RGE for A3 is proportional to A3 itself. Consequently the 
structure of the superpotential in eq. (1.6) is scale invariant. The same argument are 
valid in the SSSM with three generations. 

To establish the conditions on the soft parameters we need first make the rotation on 
the superfields that eliminates fx^ from the superpotential [31, 92] 

^0 — — ^0 -\ ^3 

L,^ — ^3 ^0 

where /^^ = /Xq -|- The sneutrino vev in this basis is just 

^3 = (1-12) 
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where = (La), a = 0,3. The bihnear term cannot be rotated away from the SUSY 
lagrangian unless that ^o/vq = fis/v^. In this case is ahgned with Va, the L violating 
bilinears terms can be rotated away from the SSM lagrangian and the neutrino acquires 
their mass only at the loop-level through the induced trilinear L violating term. However 
in general 0. The misalignment can be quantified by means of an angle ^ defined as 
[88, 89, 96, 99, 100], 

sm^ = — = (l-l>j) 

Vd fJ'Vd 

where v"^ = Vq + fg. Using the minimization equations sin^ can be written in terms 
Am2 = ml^ - ml^ and AB = - Bq [31, 37, 101] as 

. . IJ^sfJ'O ( foAm^ VuIJ,Ab\ 

sm^ = 5- TT^ ?r\ (1-14) 

/i^ \vd m'lo Vd m'lo j 

where m'~o is the tau sneutrino mass in the MSSM. Consequently the necessary conditions 
in the soft terms are Am^ = and Ai? = 0. It is worth to stress that they correspond 
to universality conditions in Supergravity (SUGRA) scenarios naturally realized only at 
the unification scale (~ 10^^ GeV). In fact, to suprcss sin^ we need to require universality 
only on lepton and Higgs soft masses. Universality will be effectively broken at the weak 
scale due to calculable renormalization effects. For definiteness and simplicity we will 
adopt this assumption throughout this work, unless otherwise stated. 

Consequently, the presence of the bilinear term L^Hu in the scalar potential in the low 
energy superpotential is unavoidable if one assumes that some L violation was produced 
in the superpotential of the theory at some large scale (such as Planck or grand unification 
scale). In particular, the effect of L^H^ in the generation of the neutrino mass can be 
combined either with the effect of L^H^ or L^Q^D^ term in the superpotential, depending 
on the basis choice. 

Therefore, one could argue that models which break explicitly L, in which //j are 
neglected, may be considered to be intrinsically incomplete, if not inconsistent. 

It is worth noting also, in contrast with spontaneous i?-parity violation ([34] and 
references therein) that doublet sneutrino vev in the bilinear model is much more loosely 
constrained because it is not subject to constraints from astrophysics [102] 

We now turn to the calculation of the neutrino mass in an basis independent way. 

As observed before, the one generation SSM can be described in various equivalent 
bases, for example 

1. one in which bilinear term and sneutrino vev are non-zero, ^ and v| 7^ [23, 
31, 62, 91, 92, 95, 103-108] 

2. one in which trilinear term and sneutrino vev are non-zero, X^^ ^ and 7^ [109] 

3. the vev-less basis in which //|^^ and Af^^ are non-zero but v^^^ — [110, 111] 
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where the L violating parameters can be expressed in terms of dimension-less basis- 
independent alignment parameters sin^ sin^' and sin^" [31, 92, 112] [X = I, II or ///) 
as follows: 

sin^ = tttll!tl^!l^_ri:L (1.15) 

,,X\X ,.X\X \II ,,I 

sine' = lh>^^^}H^ ^ >^ ^ _IH ^^^^^ 

_„,x\x,.,x\x „,j \in 

smC = ' ' t ' ' = !!l = -^ (1.17) 
Vdhb Vd hb 

where 

hb=]/xf + Xf l^=\/^if + l^f Vd=^/vf + vf, X^I,II,orIII 

(1.18) 

Note that, in the notation of eqs. (1.15)-(1.17), the parameters 113 and /iq appearing 
in eq. (1.8) should bear the superscript I. 

Of these parameters only two are independent because they satisfy 

sin = cos C sin ^ - sin ^ cos ^ (1.19) 

sin^ is easily generalized to three-generation case [88, 89, 96, 99, 100], 

.^^2 , _ 1 ^a,/3 (W - ^lpVaf _ 1 ^a,pAlp (,.2 = ^,.,. ,,2 = V 7; 

f'lli S~n 92 ~r> 92' VA*— ^af^aH-a: l^d — ^ai^a^a) 

2 /x^-uj 2 

(1.20) 

and q; = 0, 1, 2, 3. 

From now on we will work in the A3 = 0-basis, unless otherwise stated. As a result 
we will omit the label / in all the parameters associated with this basis. Note we can use 
hb instead Aq. One of the advantages in working in this basis is that the RGB's evolution 
does not induce the trilinear L violating terms neither in the superpotential nor in the 
scalar potential [92]. 

We define also the basis independent parameter 

tan/3' = ^ (1-21) 

Vq ^ v^ 

It makes sense in the A3 = 0-basis where the usual MSSM relation 

nib = ^^61^0 (1-22) 
v2 

to introduce the following notation in spherical coordinates for the vacuum expectation 
values (vev): 

vo = V sin 9 cos /3 
Vu — V sin 9 sin /3 

V3 — vcos9 (1-23) 
which preserves the standard definition tan/? = Vu/vq. 
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1.2.3 Origin of Neutrino Masses 

In this model the presence of bihnear L violating term in both the superpotential and the 

soft potential induces a mass for the z/3 at the tree level [16, 17, 22]. In order to study the 
z/3 mass it is convenient to have an analytical expression for m^^g in this limit. The tree 
level 1/3 mass may be expressed in an basis independent way as [31, 88, 92, 96, 99] 



{g^M^ + g''M2)iJ 



m^3 ~ 



2 ■ 2 /- 
Vd sm ^ 



;i.24) 



in terms of basis-independent parameters /i, and sin^ defined in eq. (1-18) and (1.15). 
The formula was given first in specific basis in for example [16, 21, 30, 113]. The second 
term in the denominator may be neglected if Af2,/U > mz, as often happens in minimal 
supergravity models with universal soft SUSY breaking terms [36, 114]. Thus one may 
obtain an estimate of the neutrino mass by keeping only the first term in the denominator. 



ni 



9 2 ■ 2 t 
Vd sm ^, 



2M 



;i.25) 



where we have used Mi — M2g''^/g'^. For sin^ f^i 1 one can easily check that m^^g could 
be as large as the direct experimental upper bound of 18 Mev [115]. However in SUGRA 
models with universality (in fact, we need to require universality only of lepton and 
Higgs soft masses) one may obtain naturally small sin,^ values, calculable from the RGE 
evolution from the unification scale down to the weak scale [37, 38, 41]. Indeed, using the 
minimization equations sin ^ was written in terms Am^ = — m^^ and AB — B^ — Bq 
in eq. (1-14). In term of the basis invariant quantities we have 



sm^ = -cos^ sm^ ( cos^— ^ + — 



m'-p Vd m 



(1.26) 



One may give a simplified approximate analytical expression for the in this model by 
solving the renormalization group equations for the soft mass parameters rn\^, "^L3> -^O) 
and i?3 in the one-step approximation. This gives [31, 37, 38, 92, 96, 101] 



sinC|Am2 



COS ^' sin ^' cos^h\ 
cos ^'sin ^' cos^hl 



ml, +MI + MI + A\ 



m'lp 



21 3 , Mu 
— ^ In 



D 



3 , Mu 
m 



Stt^ rrit 



(1.27) 



and 



sin C I AS ~ cos ^' sin ^' tan P'h^ 



D 



m~p 



3 - Mu 
m 



87r^ rrit 



;i.28) 



where we have denoted by the symbols sin^l^m^ and sin ^| as the two terms contributing 
to sin^ in eq. (1.26). In section 5.1.2 we compare these formulas with the full numerical 
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basis-independent 


Basis I: A| = 


Basis II: /i^^ = 


Basis III: ^ = 


sin^ 


sin ^' 


^3 


4 


^3 


^3 


^3 




(a) 




10-4 


10-2 


10-4 


1 


10-2 


1 


10-3 


10-4 


10-3 


10-4 


10-3 


10-2 


(b) 


10-5 


10-4 


10-5 


10-4 


10-^ 


10-2 


10-4 


10-4 


20-7 10-4 


10-3 


10-^ 10-5 


10-3 


10-2 



t in GeV 



Table 1.1: Estimated magnitude of L violating parameters required for a 1/3 mass in the eV 
range in the three bases defined before for (a) SUGRA (b) Horizontal symmetry. 

calculation and explore non-universality effects. Using these expressions and assuming 
no strong cancellation between these terms one finds that the minimum neutrino mass is 
controlled by the sin^jAm^- As a result one finds, [31, 37, 38, 92, 96] 




(1.29) 



The above approximate analytical form of the 1/3 mass is useful, as we will see later 
(e.g. eq. (5.11)) in order to display explicitly the degree of correlation between the L 
violating decays, such as ti br, with the 2/3 mass. 

The minimum value for sin C,'hb is determined by the value sin ^' and that of tan j3. 
For sin^' ~ 1 and relatively small tan/3 so that ht is perturbative, one has 

m^3 > lOKeV (1.30) 

for M2 ~ 1 TeV. In order to get smaller 2/3 masses one needs to suppress sin2 ^' additionally, 
for example to reach one electron-volt the required L violating parameters are given in 
Table 1.1. These order-of-magnitude estimates are given in terms of the basis-independent 
angles ^ and and in the relevant parameters for the three bases defined before. 

Note that whenever the parameter has two values, the first correspond to tan f3 — 2 
(the lower perturbativity limit) and the second to tan/3 = 35. In Table 1.1, sin^ was 
estimated from eq. (1.25) and sin,^' from cq. (1.29). 

Note also that the RGE-induced suppression depends basically in the hf factor in 
eq. (1.27) which is ~ 10-3 (~ 1) for small (large) tan/3. As a result the bigger the value 
of tan/3, the smaller sin^' will have to be for a fixed z/3 mass. The L violating parameters 
in the several bases were estimated from cqs. (1.15), (1.16) and (1.19). 

In eq. (1.29) we have neglected AB contribution with respect to the one coming from 
Am2. It is possible, however, that the AB term may be sizeable. In the large AB case 
then it may cancel the Arn^ contribution in sin^, leading to an additionally suppressed 
neutrino mass. As we will see, however, in SUGRA models with universal soft terms at 
the unification scale we do not need any substantial cancellation in order to obtain i/3 
masses below the electron-volt scale. 

In Horizontal models to be discussed in section 2.3.1, it is possible make a prediction 
for sin^. Instead of eq. (1.27) that is roughly of order hln^/ i^q, the Horizontal models 
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predict, see eq. (2.23) and [50, 89, 100, 116] 

sine~— (1.31) 

corresponding to no cancellation at all in the two terms contributing to sin^ in the first 
equality of eq. (1.15). Consequently, for a fixed neutrino mass, fii is required to be much 
lower in horizontal models than in the SUGRA case. In Table 1.1 we compare all the 
parameters in the various basis for the two cases. This is the reason why the limits on 
L violating parameters are stronger in horizontal models [50, 89, 100]. And as a result, 
the magnitude of L violating processes is correspondingly suppressed in these horizontal 
models. In this way the strength of L violation effects at colliders may give light on the 
misalignment origin of the SSSM. 

In summary for one generation, one consistent SUSY model with minimum field con- 
tent based only in gauge principles can be constructed. This model give rise to a tree 
level neutrino mass through AL — 1 violation and in general satisfies the bounds on L 
violation and as well as conserves B number. 



Chapter 2 

Supersymmetric Standard Models 



The field content of tlie Standard Model (SM) together with the requirement of Gsm = 
SU{3)c X SU{2)l X U{1)y gauge invariance implies that the most general Lagrangian is 
characterized by additional accidental U{1) symmetries implying Baryon {B) and Lepton 
flavor number (Lj , i = e , jj, ,t) conservation at the renormalizable level. When the SM is 
supersymmetrized, this nice feature is lost. The introduction of the superpartners allows 
for several new Lorentz invariant couplings. The most general superpotential respecting 
the gauge symmetries and with minimum field content reads 

W, =eab [h^Q-U.H', + KjoQl^A + ^mKL]Ek] (2.1a) 

W2^-eabl^oLlHl (2.1b) 

Ws^-Sam^H'^ (2.1c) 

^Sab [XijkL^L'jEk + Kj^L^Q'fik] (2.1d) 

W, =X''U,D,D, (2.1e) 



We ^e^^^SabScd^ QTQfQtLf + • • • (2.1f) 



where i,j,k,l = 1,2,3 are generation indices, a,b,c,d =1,2 are SU{2) indices, fiji^jP 
are SU{3) indices, and £ is a completely antisymmetric matrix, with £12 = 1. For later 
use, it is convenient to generalize the L violating superrenormalizable term to include the 
usual //-term: /laLaHu, where a — 0, . . . ,3, Lq = H^, A^^g = ^ojk = ^b- The 

symbol "hat" over each letter indicates a superfield, with Qi, Li, Lq, and Hu being SU(2) 
doublets with hypercharges |, —1, —1, and 1 respectively, and U, D, and E being SU{2) 
singlets with hypercharges — |, |, and 2 respectively. The couplings hu, Aq and Aq are 
3x3 Yukawa matrices, and ji^ and iii are parameters with units of mass. 

We call this model Supersymmetric Standard Model (SSM) but we will see below at 
least one initial description before "Supersymmetric" is necessary in order to construct 
viable phenomenological models. We will then suggest the Superrenormalizable Super- 
symmetric Standard Model as the minimal realization of SSM at low energies when the 
neutrino anomalies are taken into account. As it stands, eq. (2.1) has potentially danger- 
ous phenomenological consequences 
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i) The simultaneous presence of the Dimension-4 (D = 4) terms in eqs. (2. Id) and 
(2.1e) that violates Baryon Number (B) and Lepton Number (L) give rise to a very 
fast proton decay. 

ii) Flavor problem. The D = 4 Yukawa couplings in eq. (2.1a) are expected to be of 
order unity, suggesting that all the fermion masses should be close to the electroweak 
breaking scale. 

iii) If the dimension-4 B (or L) violating terms are not absent from the superpotential, 
the trilincar couplings Xijk , A^^-^ (or A'/^^,) are also expected to be of order unity, 
implying unsuppressed L (or B) violating processes. 

iv) The Dimension-5 {D = 5) B/L non-renormalizable violating couplings are also ex- 
pected to be order unity, implying a too fast proton decay. 

v) n problem. The superrenormalizable parameters Ha are gauge and supersymmetric 
invariant, and thus their natural value is expected to be much larger than the elec- 
troweak and supersymmetry breaking scales. A large value of fiQ would result in too 
large Higgsino mixing term (this is the supersymmetric n problem) . The parameters 
/Ii are required to be further suppressed by the smallness of the neutrino Masses 

vi) Large neutrino masses: the presence of general D = 4, L violating terms expected to 
be order unity imply a potentially large tree level neutrino mass. 

All these puzzles strongly indicate that SUSY models should be restricted by some addi- 
tional symmetry other than SUSY and SU (3) x SU (2) xU{l) gauge symmetry. In general 
we label such symmetry as X symmetry. 

The lack of explanation for the order unity couplings in eq. (2.1) is a problem common 
both SM and SSM. But we know SM Yukawa couplings are of order hij ~ mirrij j . 
In this way and following [11] we call natural conditions on the dimension-4 and 5 L 
and B violating couplings to require them to be order Ajjfc ~ ^J'mimjmk/v■^ and k'^^^i ~ 

^J mirrijmkmi / respectively, so that Xijk ^ 10~^ and n'ijk^ ^ 10~^. In this case the 
problems in iii) [11, 41] and vi) [96] disappear. However, proton stability forces A' and A" 
to be much smaller [24, 117] 




k 7^ 1 



(2.2) 



and [11] 



KjkKik < 10 



1-24 



J = 1,2 



(2.3) 



When higher generations are involved, weaker constraints apply [37] 



KjkKmn ^ 10 



(2.4) 



Concerning problem iv), SUSY models are sensitive to flavor physics through B vio- 
lation suppressed by Planck scale. For instance, the operator {}-/ Mp){QiQ]){Q2Li) gives 
a proton life time shorter than the experimental bound by about 14 orders of magnitude 
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[118]. An example of a contributing Feynman diagrams shown in Fig. 2. In general the 
present upper bound on the proton decay rate can be translated on k'-^-^; in the following 
way [119] 



1121 



< 10-2^ GeV 



MsusyM 



(2.5) 



If, for example, we take M ~ Mp/VStt = 2.4 x 10^^ GeV and Msusy ~ 10^ GeV, then 
eq. (2.5) become 



<121 < 



(2.6) 



P 




d 



u 



Figure 2.1: Possible contribution to proton decay, involving n'^^f^i 

Thus the simplest supersymmetric extension of the SM is excluded: an extra X sym- 
metry is required to protect the proton. Moreover stability of proton allows either L or B 
violation, because the nucleon decays normally requires not only baryon number violation 
but also lepton number violation. This implies at least the absence of either dimension-4 
L OT B violating terms and the absence of the dimension-5 B/L violating terms. The 
most popular example of X symmetry is -R-parity which forbids all the dimension < 4 B 
and L violating terms. Consequently it avoids the most dangerous constraints in eq. (2.2). 
However it was promptly realized [117] that such a symmetry emerging for example from 
minimal SU{5) SUSY [6], is in conflict with proton decay because of the limit in eq. (2.6). 
In fact, the minimal SU{5) SUSY realization of Rp [6] is nearly excluded by a combi- 
nation of the as{mz) prediction and the proton life time, and will be severely tested at 
Superkamiokande and LEP2 [120]. Since the Planck scale operators give proton decay 
rates which are too large, one needs another suppression mechanism. 

In this way the effect of the dimension-5 B/L violating operators cannot be neglected 
neither in strings nor Grand Uniflcation Theories (GUT), even after assuming that they 
have natural couplings [11]. This suggests that B and L conservation in SUSY models 
may not be a consequence of R parity, but of a different X symmetry. 

2.1 Possible X symmetries 

There is a wide variety of possibilities for the X symmetry (see [4] and references therein). 
One key ingredient of a such symmetry is to have a justified origin. The X symmetries 
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have their origin in one of the three categories of symmetries which occur in field theory 
models of particle physics: space time symmetries, gauge symmetries or flavor symmetries. 
The X symmetry most frequently used is i?-parity. Since it acts in the anti-commuting 
coordinate of superspace it can be viewed as a superspace analogue of the familiar discrete 
space-time symmetries. However the experience with space time symmetries such as P 
and CP, that in the real world are broken, suggests that broken R parity models are a 
likely possibility. 

Assuming that the low energy discrete symmetries come somehow from a larger gauge 
symmetry gives a rationale for the very existence of discrete symmetries. Otherwise, the 
existence of discrete symmetries seem unmotivated from a fundamental point of view. 

We start considering the possibility for X to be a discrete subgroup of an enlarged 
anomaly free gauge symmetry. Such a symmetry will be an anomaly free discrete sym- 
metry that could arise naturally in string theories [7-11]. This kind of symmetries include 
all previously discussed symmetries used in this context such as matter parities [121] 

2.2 Gauge Discrete symmetries 

In = 1 supersymmetric theories there are two types Abelian internal symmetries: 
ordinary symmetries and /^-symmetries. The flrst commutes with the SUSY generator 
and the second does not because the superspace grassman variable 6 is charged under the 
i?-symmetry . They lead to two types of discrete symmetries. They are (a) discrete Zjv 
symmetries [7, 8, 10, 11, 13, 40, 119, 122] or (b) discrete R symmetries [10, 41, 123-125] 
We will study in detail only the (case (a)) because the analysis for the other is 
rather similar[10. 123]. Specifically we consider a Zjv symmetry under which each chiral 
superfield transforms as 



where the are additive Z^r charges. An operator is allowed only if the sum of its 
charges is (mod A"). We will assume that the Zjq charges are not family dependent. 
A global discrete symmetry is not protected against violation by Planck-scalc and other 
non-perturbative effects (see [10]). Moreover the origin of discrete symmetries is arbitrary 
unless they are gauge discrete symmetries (see [10] and references therein). Furthermore 
it was pointed in ref. [9] that discrete gauge symmetries are restricted by certain anomaly 
cancellation conditions. Thus many candidate gauge discrete symmetries may be ruled out 
on the basis of these conditions. One way [12, 40] to obtain a gauged discrete symmetry 
is to break a gauged U{1) symmetry with an order parameter (vev) whose charge is 
normalized to g, where the smallest non-zero f/(l) charge assignment in the theory is 
q [12, 13, 126]. For a complete discussion on the possible origin of these gauge discrete 
symmetries see [10]. 

An analysis of the anomaly cancellation of discrete symmetries was given in [7, 9-11, 
123, 127]. We will follow the notation of [11] and write down the charges of the superfield 
as a vector of the form 



(2.7) 



a = {aq, a^, a^, ae, aH^, o^hJ ■ 



(2.8) 
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The presence of the Yukawa couphngs of eq. (2.1) and invariance under hypercharge, 
reduce the number of independent couphngs to just three [11]. Thus we can choose a 
convenient basis in which the charge of any field is given in terms of three integers {m,n, 
and p) 

(0,-1,1,0,1,-1,1), (2.9a) 
aA = (0,0,-1,-1,0,1,0), (2.9b) 
= (0,0,0,-1,1,0,0). (2.9c) 

So the total charge can be written as 

a — man + naA + potL — (0, — m, m — n, —n — p,m + p, —m + n, m) (2-10) 

The elements of the basis in eq. (2.9) can be considered as the three independent Zjv 
generators discussed in refs. [7, 10] 

i?Ar = exp(i27rQ;ij/A'") , Ajv = exp(i27rQ;A/-/V) , — e:K.p{i27raL/N), (2.11) 

In terms of these generators one can write the Flavor-independent discrete symmetry as 
[7, 10] 

= i?™ X A% X LI (2.12) 
The conditions to have the various terms in the superpotential of eq. (2.1) are 



W2 n = (modA^) (2.13a) 

W3 => m-n-p^O (modA^) (2.13b) 

W4 m-2n-p = (mod A^) (2.13c) 

W5 m-2n = (modA^) (2.13d) 

We -n-p^O (modA^) (2.13e) 



where we have corrected a misprint in the sign of eq. (2.13e) in [11], and the Wi were 
defined in eq. (2.1). They correspond to terms with couplings /Iq, /li, X'^j^ (and Xijk), Xll^ki 
and n',ijf^i respectively. 

If the term is allowed in the superpotential (n = modA^), the possible models 
are quite restricted. We show them in the upper part of Table 2.1. We also show there 
the equivalent matter parity names [121]. The generalized Z^v-type matter parity is also 
given according to the definition given in eq. (2.12). Finally we display the conditions on 
n, m and p that need be satisfied in order to obtain the corresponding set of operators. All 
these conditions are mod A^. The column for A'^^-^, actually includes also the Xij^. coupling. 

Note that all models in n = case forbid proton decay avoiding the simultaneous pres- 
ence of renormalizable L and B violation. Moreover, the anomaly cancellation conditions 
impose restrictions on A^, m, n, p. The corresponding possibilities are the symmetries 
[10]: R^Ll (MSSM, n = 0, m = 1, p = 2 in Table 2.1); i?2, (i?p-SSM)i; L3 (L^-SSM); 
R^L^ (5p-SSM); and RlL^ (MSSM) [128] 

^The i?3 symmetry does essentially the same job as the standard Rp, at least with respect to the 
operators displayed in Table 2.1 
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^ This is a simplifying notation. The real condition is p — m = (mod AT) . Consider for 
example: e.g, iV = 3, n = —6, m = —2, p = 4. 



* Additional condition: m^2n 
*Additional condition: p 7^ — n 
•Additional condition: m 7^ 2n 

Table 2.1: Possible models arising for discrete symmetries. All the relations in n, m, p 
shown in the table are also mod N. Note that the n = case is very restrictive in the 

type of possible models. For the case of n 7^ (modA^) the simultaneous presence of 
superrenormalizable and dimension-4 L violating operators is not allowed. 



2 Supersymmetric Standard Models 



21 



As the fiQ is likely to have a different origin than the Yuk;awas in Wi we also will 
consider the possibility that the fiQ can be absent from the high energy superpotential. 
This corresponds to the extra lines with n 7^ in Table 2.1. In this case it will be 
necessary to use a discrete version of the Green-Schwartz mechanism (GS) [129] , and the 
number of potentially anomaly free solutions substantially increases. Therefore, this kind 
of symmetries may be well suited for the solution of the fi problem. In fact, once mass 
parameters are introduced, the term fiaLaHu (rather than only ^^LdHy) is the most 
natural choice. Consequently the SSM model with minimal set of couplings in this case, 
is likely to include general superrenormalizable terms but not dimension-4 ones. We call 
it Superrenormalizable Supersymmetric Standard Model (SSSM). The main difference of 
the SSSM with the general i?p-SSM is that in the former case the induced dimension-4 L 
violating couplings are naturally suppressed and automatically satisfy the experimental 
constraints on L violation. Examples of anomaly free (through GS) discrete symmetries 
in this case are [10]: Z2 symmetries A2L2 {n — 1, m — 2, p — 1) corresponding to A'}^L% 
case (m = 2n, p = n, 2n = 0) in Table 2.1. and R2A2L2 {R^^^^A^'U' in Table 2.1). 
symmetries with > 3 may be made anomaly free now. Examples are [10]: A^^L^ that 
correspond to A^}^L~^ case in the table (SSSM plus QQQL); and A^L^ that correspond 
to R^^A^L"^ case in Table 2.1 (SSSM) 

The possibilities for anomaly free gauge discrete symmetries stabilizing the proton 
further increases if one allows for gauged it!-symmetries [10] 

The resulting models compatible with proton decay are hsted in Table 2.2. 



Dim. < 5 L OT B violating terms 


Example 


Symmetry 
Name : type 


Origin 


A 


[128] 


MSSM : RlLs 


Strings 




[8, 52] 


SSSM : A^Ll 


Strings 




[119] 


LpSSM : L3 


Flavor 


\jkLiLjEk + ^iLiH^ 


[130] 




GUT 




[43] 




GUT 


KjkLiLjEk + X^ji^LiQjDk + fXiLiHu 


[7] 


BpSSM : R3L3 


Flavor 



Table 2.2: Possibilities of low energy Lepton or Baryon number violating superpotential 
renormalizable terms where the proton should be enough stable. For each case one specific 
model is cited with their correspondent discrete symmetry if possible. 

The triangle in Table 2.2 is just to illustrate that when we roll down the "hill" the 
number of parameters in the superpotential increases: 0,3,9,12,30 and 39 respectively. 
This is the main reason why most of the studies were initially performed in the context 
of the MSSM (really in the context of i?p-SSM). When the evidence of neutrino masses 
and mixings is taken into account clearly the minimum model which emerges is just the 
Superrenormalizable SSM. We now review all the possibilities with emphasis in the SSSM 
alternative. 
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2.2 Gauge Discrete symmetries 



We define the superpotential with minimum content of couplings as the one without 
any dimension < 5 L/B violating operators [10, 128] 

W^MSSM = W, + W2 (2.14) 

On the other hand, i?-parity is equivalent to the R2 anomaly free gauge discrete symmetry 
{m — 1, n — p — 0) 

i?2 = -R3 =^ {o!Q, ««, (Xd, OLL-, ae, ftij^, anj = (0, -1, 1, 0, 1, -1, 1) . (2.15) 
In the i?2 case the high energy superpotential contains [128] 

Wr^ = Wmssm + '^Q^QJQkLl + KLHuS + MsSS (2.16) 

The latter allows both proton decay and neutrino masses with Rp conservation. For R3 we 
should have XsSSS instead of the last term. The introduced massive Majorana superfield 
S with i?2 charge as = —1, is necessary to cancel the mixed gravitational anomaly [10] 
and also leads to a AL = 2 Majorana neutrino mass. 

The presence of the D = 5 term may give to rise to fast proton decay [4, 7, 10-12, 117, 
119, 123]. 

2.2.1 Superrenormalizable Supersymmetric Standard Model (Bilinear Model) 

From now on we focus on the i?p-SSM superpotential defined in eq. (2.16). There are 
two main approaches to break Rp in the bilinear way, explicitly and spontaneously. See 
Fig. 2.2. 

In the spontaneous approach, the breaking of Rp occurs by minimizing the scalar po- 
tential and the most relevant parameter is the vacuum expectation value of the sneutrino 
which breaks lepton number [14, 15, 17, 27, 52-54, 56, 58-60, 131, 132]. 

For the ungauged L case (see Fig. 2.2) there is a Goldstone boson (Majoron) so 
that the theory must contain a singlet right handed neutrino superfield [52, 53] with a 
vcv which typically lies at the weak scale. For the gauged case the Goldstone boson is 
absorbed by a new gauge boson [51, 133]. Either way one arrives to a superrenormalizable 
supersymmetric model with a relatively small /li, as indicated by present anomalies in 
physics. 

Here I will concentrate mainly on the explicit Rp violating case in which the super- 
renormalizable nature of the Rp violation mechanism must come from new physics near 
the Planck scale (see section 2.3). Note that in an explicit Rp violating model (like 
the SSSM) one can have also a sneutrino vev. Conversely a spontaneous model may be 
accompanied by small explicit Rp breaking terms. 

The resulting renormalizable superpotential which emerges is of the form 



W^SSSM = Sab 



KQWiK + h%Q\D^Ll + h%L\E^Ll - pt^LlH', , (2.17) 
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the soft superpotential is expected to satisfy the same symmetry of the superpotential 
and therefore we have 

V^oft = MfQTQ"; + MfU*Uj + MfD*Dj + Mf^L''^*^ + MfE*Ej 
- [iMgAsAs + \M\2\2 + |MiAiAi + h.c] 



-E-^L^ifg , (2.18) 

If the non-diagonal entries of M^^ and Sq/j turn to be neghgible we can simply define: 
mn, = Mf , Ml, = Bo = Boo and B^ = B^. 

In SSSM the appearance of misalignment between the Va = {La) and fia, param- 
eterized by sin^ in eq. (1.20), is unavoidable with the possible exception of one point 
characterized by very specific relations between some of the soft parameters. In the one 
generation case they were simply Am^ = and AB = in eq. (1.26). In the three- 
generation case these conditions are [16, 89, 92, 96, 100]. 

(A) /Iq, is an eigenvector of e.g, m|^^/iQ, = mlfia 

(B) Ba/sfi/s/fia = B 6a/3 for all a. 

and are naturally realized by the universal conditions at the unification scale. However, 
once sin^ 7^ 0, radiatively induced for example, the term proportional to L^Hu cannot 
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be rotated away of the SSSM lagrangian. They give rise to a plethora of new effects 
without counterpart in a model with only dimension-4 L violating terms (absent if only 
dimension-4 L violating terms were considered) like the emerging a tree level neutrino 
mass, additional contribution to radiative neutrino masses, new decay modes for superpar- 
ticlcs, sncutrino splitting, and so on. See section 2.4 for a review on the phenomenology 
of SSSM. Unless otherwise stated, we will assume universality of at least lepton and Higgs 
soft masses at the unification scale in the SSSM. In this way we may obtain naturally 
small tree level neutrino mass, induced by the RGE running from the unification scale 
down to the weak scale as explained in section 1.2.3. 

We consider the model defined in eq. (2.17) and eq. (2.18) as our reference Superrenor- 
malizable Supersymmetric Standard Model. 

Note that L violation can be rotated away from the superrenormalizable L violating 
terms (but not from bilinear soft terms) into L violating dimension-4 terms 

= ^h%L,QjD, + f^hi%LjE, (2.19) 
A* A* 

where n"^ = Hq + n]. The new L violating couplings are not arbitrary and lead the same 
phenomenology of the general Bp-SSM (see [87] for a detailed discussion on this). Note 
also that the rotation leading to (2.19) is not scale invariant and that W^^^ will regenerate 
L violation in superrenormalizable terms via renormalization at scale A'. So that a model 
without bilinear L violating terms is incomplete and inconsistent. 

Before closing this section let us comment that another way to generate the SSSM is 
provided by the first row with n 7^ in Table 2.1 [40, 41, 50, 123]. We will study one 
specific example [50] in the next section. 

2.3 Alternative SSSM origin 

Models about the origin of SSSM (with explicit fii terms) have been constructed in the 
context of Strings [8, 42]; GUTs [16, 37, 38, 48, 49, 130, 134-136]; i?-symmetries [41, 47, 
137]; Peccei-Quinn symmetries [39, 45, 46, 138]; and fiavor symmetries [40, 44, 50]. Next 
we will consider the latter possibility. 

2.3.1 Flavor symmetries 

The emergence of string theories as a universal theory encompassing all known fundamen- 
tal interactions including gravity provides a framework which allows to relate features of 
the effective low energy theory which seemed otherwise uncorrelated. Of special interest 
for the problems that we are discussing here is the presence of non-renormalizable inter- 
actions, coming from a large number of horizontal gauge symmetries, especially Abelian, 
which are spontaneously broken at scales that may vary between the electroweak scale 
and the Planck scale. 

All these properties lied as to reconsider the original idea of Froggatt and Nielsen (FN) 
[139]. The approach originally suggested by FN to solve the fiavor problem and account 
for the fermion mass hierarchy, turns out to be quite powerful in the context of the SSM to 
solve also the /i problem. FN postulated an horizontal U{1)h symmetry that forbids most 
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of the fermion Yukawa couplings. The symmetry is spontaneously broken by the vacuum 
expectation value (vev) of a SM singlet field x ^-nd a small parameter of the order of the 
Cabibbo angle 6 = (x)/M ~ 0.22 (where M is some large mass scale) is introduced. The 
breaking of the symmetry induces a set of effective operators coupling the SM fermions 
to the electroweak Higgs fields, which involve enough powers of 6 to ensure an overall 
vanishing horizontal charge. Then the observed hierarchy of fermion masses results from 
the dimensional hierarchy among the various higher order operators. When the FN idea 
is implemented within the SSM, it is often assumed that the breaking of the horizontal 
symmetry is triggered by a single vev, for example the vev of the scalar component of a 
chiral supermultiplet x with horizontal charge H{x) = — 1 ■ 

More recently it has been realized that the FN mechanism can play a crucial role 
also in keeping under control the trilinear B and L violating terms in (2.1) without the 
need of introducing an ad hoc R-parity quantum number [12, 13, 25, 40, 44, 50, 116, 126, 
140, 141]. There are three ways to naturally suppress the B and L violating couplings in 
models with horizontal symmetries without require too large horizontal charges for the 
SSM fields 

(a) Since H{x) — —1 , then because the superpotential is holomorphic all the operators 
carrying a negative charge are forbidden in the supersymmetric limit. [13, 44, 50, 116, 
126, 141] 

(b) If the operator carry a fractional charge only can appear at very high dimensional- 
ity [40, 44, 50, 116, 126] 

(c) If the horizontal symmetry is not completely broken: a residual symmetry remain 
that forbids the dangerous operators [7, 12, 13, 40, 44]. In fact a simple way to obtain 
either Bp [13] or Lp [12] is by the spontaneous breaking U{1) which arises if 
the field which breaks U{1) has a charge N normalized to the smallest charge of the 
theory 

The mechanism (a) is very related to the solution of the /iQ-problem. If under U{1)h 
the superrenormalizablc terms LqH,^ has a charge no < , a /io term can only arise from the 
(non-holomorphic) Kahler potential, suppressed with respect the supersymmetry breaking 
scale 777.3/2 as [142] 

/io ^ ms/2 . (2.20) 

A too large suppression (|77o| > 1) would result in unacceptably light Higgsinos, so that 
in practice on phenomenological grounds no = —1 is by far the preferred value. 

For example in [141] it was argued that under a set of mild phenomenological assump- 
tions about the size of neutrino mixings a non-anomalous U{1)h symmetry together with 
the holomorphy conditions implies the vanishing of all the superpotential B and L violat- 
ing couplings. A systematic analysis on the restrictions on trilinear L violating couplings 
in the framework of U{1)h horizontal symmetries was also recently presented in [140]. 

In this work we argue that if the yUo problem is solved by the horizontal symmetry in 
the way outlined above, and if the additional bilinear terms /ij are also generated from 
the Kahler potential and satisfy the requirement of inducing a neutrino mass below the 
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eV scale, as indicated by data on atmospheric neutrinos [26, 143, 144], then in the basis 
where the horizontal charges are well defined, all the trilinear L violating couplings are 
automatically absent. This hints at a self-consistent theoretical framework in which L 
is violated only by bilinear terms that induce a tree level neutrino mass in the range 
suggested by the atmospheric neutrino anomaly, L and B violating processes are strongly 
suppressed, and the radiative contributions to neutrino masses are safely small so that 
^loop ^ 10^^ eV, which barely allows for the LOW or quasi-vacuum solutions to the solar 
neutrino problem [26, 145, 146] with tan/3 in the range ~ 10-40. 

In the following we will denote a field and its horizontal charge with the same symbol, 
e.g. H{ei) = Ci for the lepton singlets, H{Qi) = Qi for the quark doublets, etc. It is 
also useful to introduce the notation fij = fi — fj to denote the difference between the 
charges of two fields. For example Lio denotes the difference between the charges of the Li 
'lepton doublet' and the Lq 'Higgs field'. On phenomenological grounds we will assume 
that the charge of the /iq term is no = —1 and we will also assume negative charges 
rii — Li + Hu < riQ for the other three bilinear terms LiHu ■ It is worth stressing that the 
theoretical constraints from the cancellation of the mixed Gsm x U{1)h anomalies hint 
at the same value no = — 1 both in the anomalous [147] and in the non-anomalous [141] 
U{1)h models. With the previous assumptions the four components of the vector in 
(2.1) read 

where coefficients of order unity multiplying each entry have been left understood. It is 
well known that if Ha and the vector of the hypercharge Y = —1/2 vevs Va = (La) are 
not aligned [16, 100] (sin^ 7^ in eq. (1.20)) the neutrinos mix with the neutralinos [16, 
22], and one neutrino mass is induced at the tree level giving by eq. (2.21) which can be 
rewritten as 



where = Micos'^Ow + M2 sin^ , Mi and M2 are the f/(l)y and SU{2)l gaug- 
ino masses, and tan/3 = {Hy) / {Hd) ■ Since nib/ {Hu) tan j3 ~ 6'^'''tan/3 (with mb{mt) ~ 
2.9 GeV [148]) in the following we will use the parameterization tan P — 9^~^ that ranges 
between 90 and 1 for x between and 3. Keeping in mind that we are always neglecting 
coefficients of order unity, we can approximate cos^/9 = (1 -\-ia.v? (3)~^ ~ ^2(3-2:)^ Taking 
also Ml ~ , 11M2/MI > sin 2/3 cos^ and 100GeV< M2 < 500 GeV we obtain from 
(2.21) 

^tree _ -{5+x) ^] %V . (2.22) 

The magnitude of the tree- level neutrino mass as a function of logg sin ^ ^ L30 for different 
values of x (which in our notations parameterizes tan/3) is illustrated in fig. 1. The grey 
bands correspond to equation (2.21) with M2 ranging between 100 GeV and 500 GeV, 
while the dashed lines correspond to the approximate expression (2.22). 

In general, the two conditions (A) and (B) in section 2.2.1 have to be satisfied to ensure 
exact Ha^Va alignment and m^^^ = . In our case the goodness of the alignment between 
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Figure 2.3: Tree-level neutrino mass dependence on log^sin^ ^ L30 for different assign- 
ments of the charge difference L30 = — Lq and for different values of tan (3. Details in 
the text. 

fia and Va is controlled by the horizontal symmetry, and in particular there is no need of 
assuming universality of the soft breaking terms to suppress w}J'^'^ to an acceptable level. 
This is because the previous two conditions are automatically satisfied in an approximate 
way up to corrections of the order ^l-^^l , where the minimum charge difference between 
Lq and the Li 'lepton' fields is responsible for the leading effects. Thus we can estimate 



Sin^ ~ 0\L^o\ — 0\ni~no\ ^ /fi _ 



(2.23) 



Confronting (2.23) with (2.22) it follows that in order to ensure that m^J'^'^ is parametrically 
suppressed below the eV scale we need 



lui — no\ > 5 + X 



1,2,3) 



Let us introduce the parameterization 

\ni - no\ - {5 + x) = 6i . 
Without loss of generality, we can also assume ni < n2 < which implies 

It is worth stressing that the parameter that controls the scaling of m}^'^'^ with respect to 
changes in the values of the horizontal charges is ~ 0.05, and thus neutrino masses 
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are much more sensitive to the horizontal symmetry than the other fermion masses that 
scale with 9. For example ^3 = —1 yields rn"^'^'^ ~ 20 eV; ^3 = yields m}^^^ ~ leV; 
for ^3 = 1 all the trilinear L violating couplings are forbidden, and at the same time 
m}'^^ ~ 5 X 10"^ eV (see Fig. 1) is in the correct range for a solution to the atmospheric 
neutrino problem [26, 143, 144]; finally, 5^ — 2 would suppress w}^^ too much to allow 
for such a solution. 

Let us now write the down-quarks and lepton Yukawa matrices as 



IT'jk 



where x — L^ + Qz + d^ — Lq + L^ + e^ consistently with our parameterization of tan/? and 
with the approximate equality between the bottom and tau masses at sufficiently high 
energies (which in particular allows for 6-r Yukawa unification). The order of magnitude 
of the trilinear L violating couphngs is then: 



A' 



ijk 



j3+<ifc3-(5+<5i) 



(2.24) 



One can show that the phenomenological information on the charged fermion mass ratios 
and quark mixing angles can be re-expressed in terms of the sets of eight charge differences 
given in Table 2.3 [40, 116, 149-151]. Consequently gives rise to eight conditions on the 
fermion charges 



model 


Ql3 Q23 di3 0?23 ^^13 U23 


model 


-^13 + ^13 L23 + 623 


MQl: 
MQ2: 


3 2 1 5 2 
-3 2 7 11 2 


MLl: 
ML2: 


5 2 
9 -2 



Table 2.3: Models of horizontal charge differences 



The phenomenological analysis leading to these sets of charge differences has been 
extensively discussed in the literature [40, 116, 149-151]. 

In [50] was shown that in the framework of models of Abelian horizontal symmetries, 
the phenomenological information on the charged fermion mass ratios and quark mixing 
angles expressed in terms of the eight horizontal charge differences in Tab. (2.3), when 
complemented with the requirement that m]J^'^ is adequately suppressed below the eV 
scale {6i > 1) hints at one self-consistent model (MQl-l-MLl) where all the A and A' 
couplings vanish. It is interesting to note that ^3 = 1 which yields m}J^^ ki Q"^ eV in the 
correct range required by the atmospheric neutrino problem is also the minimum value 
that ensures A = , A' = and, as we will below. A" = . 

Concerning MQ2 either the neutrino masses are uninterestingly small there, or the A' 
conflicts with existing experimental limits. MQ2 is also excluded by the requirement that 
the A" couplings vanish 
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In ML2, once we set L23 = to allow for maximal v^-Vt mixing, the lepton mass ratios 
can be correctly reproduced only if L12 > 4 , which would again exclude the possibility of 
explaining the solar neutrinos deficit through Ue-i'jj, oscillations. 

In conclusion, our analysis results in the set of fields charge differences and of Ua — 
La + Hu charge sums displayed in Tab. 2.4 



Ql3 Q23 di3 d23 Ui3 U23 


-^13 L23 ei3 623 


rii no 


3 2 1 5 2 


5 2 


-7-x -1 



Table 2.4: Final charge differences 



If we further use the analysis of the neutrino loop effects we would have in addition 
the value x — 1 (corresponding to tan/3 ~ 10-40). Of course, x = implies that the 
value of tan /3 is very large (~ 60) and therefore this case is phenomenological disfavored 
[152, 153]. Finally x = 2 (tan/5 ~ 5) would yield a too large suppression to the loop 
neutrino mass to be interesting for the solar neutrinos. 

Models based on a single U{1)h Abclian factor arc completely specified in terms of the 
horizontal charges of the SM fields. There are five charges for each fermion family plus 
two charges for the Higgs doublets, for a total of 17 charges that a priori can be considered 
as free parameters (the charge of the U{1)h breaking parameter 9 is just a normalization 
factor). The 17 horizontal charges are constrained by eleven phenomenological conditions 
corresponding to the eight constraints in Table 2.3, the rrit value, the condition rrib ~ rrir 
at high energy and the preferred value for the solution of the Hq problem hq = —1. Also 
by two theoretical conditions from anomaly cancellation through the GS mechanism. 
This leaves us with four free parameters, and we can chose them to be the charges Ui 
{i = 1, 2, 3) of the bilinear terms Hi, and x = Q3 + + Lq that fixes the value of tan (3 . 
The expressions of the horizontal charges for all the SM fields as a function of these 
four parameters is given in the Appendix A. Then the self-consistent solution given in 
Table 2.4 is used to find the 17 individual charges presented in Table A.l. In the present 
case the Kahler contributions to the trilinear couplings are of order m3/2/Mp and powers 
of 6 depending on the horizontal charges of the various dimension-4 terms. Moreover, 
for this set of individual charges the horizontal charges of A" are in addition fractional 
and consequently zero. In summary the resulting Dimension 4 and 5 L and B violating 
couplings in the basis where the horizontal charges are well defined are 

/i. = m3/2^^ ~ (m3/2/Mp)^l''^+'^^+''=l, Xijk ~ (m3/2/iWp)^l^'+^^+^'=l, 

Kjk — 0> l^'ijkl — 0- 

(2.25) 

This clearly define one SSSM. 

In summary we have obtained one SSSM originating from an anomalous horizontal 
symmetry where the anomalies are canceled through a GS mechanism. We have assumed 
that all the superrenormalizable terms coupling the up- type Higgs doublet with the four 
hypercharge —1/2 doublets carry negative horizontal charges, and hence are forbidden by 
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holomorphy. We have constrained the value of these charges by several theoretical and 
phenomenological requirements, such as having an acceptable Higgsino mass {fi problem) 
and neutrino masses suppressed below the electron-volt scale, as suggested by present 
neutrino data. We have found that under these conditions all the trilinear L violating 
superpotential couplings vanish, yielding a SSSM which is defined by the charge differences 
in Tab. (2.4), where lepton number is mildly violated only by small bilinear terms. The 
model allows for neutrino masses in the correct ranges suggested by the atmospheric 
neutrino problem and by the LOW and quasi-vacuum solutions to the solar neutrino 
problem. However, no precise theoretical information can be obtained about the neutrino 
mixing angles except for the fact that, unlike the quark mixings, there is no parametric 
suppression of their values and thus they can be naturally large. Note that this model 
solves all the problems mentioned at the beginning of this chapter. 

A similar mechanism but without the solution to the /j, problem was implemented 
in [44] . In the example considered there the dimension-4 L and B violating operators are 
absent by holomorphy before the SUSY breaking. The dimension-4 L and B violating 
terms are generated from the Kahler potential and one get A'A" ~ 10^26^ However the 
strength of the superrenormalizable term 9^~^^^ is too small to be relevant in neutrino 
physics. The resulting individual charges for pa m3/26'l~^l are shown in Table A. 2 . The 
dimension-4 L and B violating terms satisfy now A'A" < 10~^^. 

2.4 Phenomenology of SSSM 

The phenomenology of the SSSM has been extensively studied in literature. In particular: 

• Neutralino-neutrino mixing: The tree level neutrino mass has been studied in [16, 

18, 19, 100, 154, 155]. Analysis of one-loop neutrino masses and/or Solution to the 
neutrino anomalies can be found in [16, 88, 89, 94, 156-165]. Numerical calculation 
with the full one-loop corrections: [34, 38]. An analytical study including the efi'ect 
of arbitrary dimension-4 L violating terms was performed in [96]. The neutrino 
spectrum in horizontal models have been estimated in [40, 44, 50, 89, 100, 126, 140, 
158] 

Neutralino decays: [16, 19, 37, 41, 87, 88, 103, 166, 167]. Tevatron signatures of 
neutralino decay in [33]. Neutralino decay processes based on the assumption that 
the atmospheric neutrino masses and mixings are mainly due to the bilinear terms 
in [64]. In [16, 47] is concluded that the dominant decay of the neutralino is into bb 
quark jets and neutrino (or btl'^ if energetically allowed). In [33, 37, 61] it is shown 
that the decay of the lightest neutralino Xi '^iU produce comparable numbers of 
muons and taus as a result of the large mixing implied by the atmospheric anomaly. 
This result is also obtained in [34, 167]. Some of the neutralino decays channels 
for the one generation case were studied in [168]. The full numerical calculation of 
the neutralino decay in the one generation case was presented in [103], while the 
full numerical calculation in the three-generation case was presented in [167]. A 
systematic study of the SSSM including approximate formula for the neutralino- 
neutrino mixing was performed in [34, 99, 104, 167, 169]. In [170] it is shown that 
neutralino mix only with the heaviest neutrino in SSSM. 
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• Neutral Higgs - sneutrino mixing. Sneutrino production and/or decays have been 
studied in [41, 84, 87] and recently in [171]. Sneutrino mass splitting have been 
analyzed in [110, 171, 172]. An basis-independent analysis of this problem was 
performed in [173]. In [97] is studied the effect of -Ls-f^u term of the low energy 
scalar potential in the generation of spontaneous CP violation. Scalar resonance 
enhancement in e'^e~ — > Z'^Z^ W^W~ and in e~^e~ — > tt due two the sneutrino- 
Higgs mixing, have been recently studied in [174]. 

• Charged Higgs - stau mixing. Charged scalar decays have been studied in [85] 

• Stop decays have been studied in [31, 70, 105, 106] 

• Rare decays: 6 — S7 [91]; Z^ decays as Z^ ^3X1 ^-iid Z^ — * r^Xi [60, 107]; 
flavor changing coupling of the W and Z [87]; top decays [62]. The double beta 
decay have been studied in [104, 169, 175]. Proton decay in bihnear models were 
studied in [170, 176, 177]. In particular in [176, 177] the limit in eq. (2.3) translate 
into 

^Ki2 < 10-'' (2.26) 

which is important in SSSM where the dimension-4 B violating couplings are gen- 
erated from the Kahler potential without further suppression, A^'^j^, < Aijkms/2/Mp 
{j = 1,2). — s> 67 and the numerical study of misalignment [98]. — > 67 and other 
rare processes [18, 19, 178]. Leptonic phenomenology and constraints [179]. 

• Others: Effects with large neutrino mass: [108]. SSSM in gauge mediated SUSY: 
[46, 64, 158, 180]. Fermion dipole moment: [181]. Magnetic moment of the neutrino 
[182]. 

• Spontaneous breaking of R-parity. In addition to the papers quoted in Section 2.2.1, 
the phenomenology of this models have been studied in [55-57, 70, 183, 184] 
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Chapter 3 

The mass spectrum of the SSSM 



Here we discuss the tree level structure of the supersymmetric masses and mixing of 
SSSM, discussing the fermion (chargino, neutralino) and scalar (squarks, sleptons and 
Higgs) mass matrices 



3.1 Chargino mass matrix 

In the SSSM the charginos mix with the charged leptons forming a set of five charged 
fermions Ff^ i = 1, ... ,5 in two component spinor notation. In a basis where ■ip^'^ — 
[— iA"*", e^, /xj, Tj^], the mass terms in the lagrangian are 



SSSM 



' Mc^ 






Mc 




r. 



+ h.c + 



where the chargino/lepton mass matrix is given by [34] 

M2 -^Qv„ 



1 



(3.1) 



where the index i expands over the columns and k labels the rows. We note that the 
chargino sector decouples from the leptonic sector in the limit jii = Vi = 0. As in the 
MSSM, the chargino mass matrix is diagonalized by two rotation matrices U and V 



U*McV ^ = diag (m^±,me,^ = diag {mp±) 



a — 1, . . . ,5 



where i = 1,2 for the charginos, and j = 1, 2, 3 for the neutrinos. 
In the one generation case, the composition of the tau is given by 



Tr 



where ip 



+T 



iX+,HI,t'^+) and V' 



-T 



i-iX- 



72 0- 
^0^ 



(3.2) 



(3.3) 



The two-component Wcyl 



spinors and are weak eigenstates, while and are the mass eigenstates. It 
follows easily from eq. (3.2) that the matrix McM^ is diagonalized by U and the matrix 
M^Mc is diagonalized by V. 
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3.2 Neutralino mass matrix 



3.2 Neutralino mass matrix 

In the basis ^'^^ = [—i^', — ^A^, i/^, u^,, z/^, z/^] the neutral fermion mass matrix Mn is 
given by [34]: 



M 



N 



Ml 





\g' Vu 


-\g'^ 





M2 


-\gvu 


\g^ 











1 / T 


1 T 
29^ 


T 






(3.4) 



where v (/Li) is the vector defined by (fJ-a)- This matrix is diagonahzed by a 7 x 7 
unitary matrix N, 

N* Mn — diag (jn^q,mu^^ — diag (mpo a — 1, . . . , 7 j (3.5) 
where i = 1, ■ ■ ■ ,4 for the neutrahnos, and j = 1, ■ ■ ■ ,3 for the neutrinos. 



3.3 The squark mass spectrum 

The up and down-type squark mass matrices of our model have already been given pre- 
viously in Ref. [95]. Here wc generalize those to the three- generation case [31]. The mass 
matrix of the up squark sector follows from the quadratic terms in the scalar potential 



quadratic 



u 



u 



R 



Ml 



+ 



(3.6) 



given by 



(3.7) 



where Aul = l{g^ - y ) K -^1 + '^1)1 and Aur = \g' {vl -vl + vfjl are the 
splitting in the squark mass spectrum produced by electro- weak symmetry breaking, and 
A^^. = A'^h^^. The eigenvalues of are 



diag{mi,i,mii2, . . . .rriu^] 



UL 



UR 



Ml 



^ UL 
UR. 



(3.8) 



This way the six weak-eigenstate fields un and UiR {i = 1,2,3) combine into six up-type 
mass eigenstate squarks Wfe as follows: Uii = ^u^Uk = V'^^luk, UiR = ^uRUk = T*j^'^Uk. 
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For completeness, we also give the mass matrix of the down squark sector. The 
quadratic scalar potential includes 



quadratic 



Ml 



dL 



(3.9) 



given by 



(3.10) 



where Adl = "K^" + W W -^1 + ^1)1, Adr = -^^9 C^i -vl + vl% and A% 
A^^h^l). The eigenvalues of are 

diag {mj^ , mj^ , . . . , m<i J = ^ ^ 



' DL 

dr} 



(3.11) 



One is left with six mass-eigcnstatc down squarks fields dk related to diL and diR fields as 
follows: diL = Tfjk = r*^idk, diR = rl;'^4 = TtRdk. 

3.4 Scalar Mass Matrices 

The electroweak symmetry is broken when the Higgs and slepton fields acquire non- 
zero vevs. These are calculated via the minimization of the effective potential or, in the 
diagrammatic method, via the tadpole equations. The full scalar potential at tree level is 



total 



E 



dW 



+ + K 



MSSM 
soft 



+ K 



SSSM 
soft 



(3.12) 



where Zi is any of the scalar fields in the superpotential in eq. (1.31), Vd are the D-terms, 
and Vf^i^ is the L violating part of eq. (1.21) 

The tree level scalar potential contains the following hnear terms 

linear = ^dPd + ^u'^u + ^1^1 + + ^3^3 > \^-^^) 

where the different are the tadpoles at tree level. They are given by 

= -1^^/3(5/37 + 57/3)/^7 + + /^/3/^/3 " |(^^ + 9''^){v (iV (i - V^) 

(3.14) 

Repeated indexes /?,7 in eq. (3.14) implies summation over 0, 1,2,3. The five tree level 
tadpoles are equal to zero at the minimum of the tree level potential. 

We assume that non-diagonal M'j^,, soft parameters are negligible. Therefore we fur- 
ther denote as M? 
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The quadratic scalar potential includes [85, 105]. 



V, 



quadratic 



< 





HI 


s± 


HI 

Kl 




p+ 















^\ H 



~4 



+ 



(3.15) 



We will denote the physical scalar bosons by Sf, the pseudo-scalar by and the 
charged bosons by Sf . Higgs bosons mix with charged sleptons and the real (imaginary) 
parts of the sneutrino mix with the scalar (pseudoscalar) Higgs bosons. The mass matrices 
in eq. (3.15), will be given in following subsections. 



3.4.1 Chcirged scaleirs 

The mass matrix of the charged scalar sector follows from the quadratic terms in the 
scalar potential in eq. (3.15). For convenience reasons we will divide this 8x8 charged 
scalar mass matrix into blocks in the following way: 



Mi. 



- 

Ml 



(3.16) 



Using the tadpole equations given in eqs. (3.14) we obtain (sum upon repeated indices): 



This matrix reduces to the usual charged Higgs mass matrix in the MSSM when we set 
Vi — Hi — 0. The slepton block is given by 



M^ 



ML 
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k 



Mlr =4= (^^oA^* - ^loVuW^ 



x/2 



'v 4 
where = A^/i^. 

The mixing between the charged Higgs sector and the stau sector is given by the 
following 6x2 block: 



4 



The mass matrix in eqs. (3.16), is diagonalized by rotation matrices which define the 
eigenvectors 

S+ = Rs±S'+ 

and the eigenvalues in the one generation case are 

diag(0,m^±,m^±,m^±) = Rs±M|±RJ± 

3.4.2 CP-Even Neutral Scalars 

The neutral CP-even scalar sector mass matrix in eq. (3.15) is given by 



where 



M 



ss 



Mis Mi,^ 
Mi,/ M?^,^ 



Vk . td 



^0 ^0 ^0 



(3.17) 



D "^0,122, u '^'^ I 



k=l 



(3.18) 



Mi., 



-BiHi - \glvuVi 
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and 



+ bW, + + i (M^,, + Ml^,) (3.19) 



where we have defined g\ = g^ + g'^- In the upper-left 2x2 block, in the limit Vi = Hi = 0, 
the reader can recognize the MSSM mass matrix corresponding to the CP-even neutral 
Higgs sector. To define the rotation matrices let us define the unrotated fields by 



^0 ~R r.R 



R 



Then the mass eigenstates are Sf given by 

S? = Rf'S'' 



with the eigenvalues diag(m|^, . . . , = M|o (•^^*') 

3.4.3 CP-Odd Neutral Scalars 

The quadratic scalar potential includes 

where the CP-odd neutral scalar mass matrix is 



Mi.p 



(3.20) 



(3.21) 



(3.22) 



where 



BollQ llQ } ^IJ'k \ 



Bq^o— + y\ Bk^k— + — 



BoiiQ 

3 



fc=l 



Vk 



(3.23) 



Note that in the limit ^ a Goldstone boson appear. 



BiHi 



and 



2 

Lik 



1-3 



k=\ 



k=\ 



(3.24) 
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The neutral pseudo-scalar mass matrices are diagonalized by the following rotation ma- 
trices, 



(3.25) 



with the eigenvalues diag(m^j , ■ ■ ■ , n^As) — (^1^^°) > where the unrotated fields 



are 



p/O ^ 



'Pd 'pI 4 4 



(3.26) 



May be convenient some times take the pseudoscalar mass eigenvalues as input, it 
is worth briefly repeating here the discussion of the pseudoscalar bosons masses. The 
pseudo-scalar mass matrix in the one generation case is given by by: 



Bol^o^ - //o/^3^ Bono 

+/^0/^3 +B3II3 



+-63/^3 



(3.27) 



As expected, this matrix has zero determinant, since the neutral Goldstone boson eaten 
by the Z is one of the corresponding states. Therefore, the masses of the two physical 
states are given by the formula: 

m2,3 = ^IVM ± ^^(IVM)' - 4(MnM22 - Mf^ + MnMga - Mfg + M22M33 - MI3). 

(3.28) 

Therefore we can easily take one these masses as input and calculate -Bq/^o from it using 

-rri^for,,!':^ + B,.[iy,fioVy,(vo + f't + ''?) + /':-!"'y,o[/'o?'»(''o + ''D] - By,i\t(vl + t|) 



Bquo 



B3 is obtained from the minimum equation for given 113 and ^3 [85]. 
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Chapter 4 
Sparticle decays 



4.1 Decays of neutralinos 

If unprotected by the ad hoc assumption of R-parity conservation the LSP will decay as 
a result of gauge boson, squark, slepton and Higgs boson exchanges. 

Although for the discussion of flavor-changing processes, such as neutrino oscillations 
involving all three generations, it is important to consider the full three-generation struc- 
ture of the model, for the following discussion of neutralino decay properties it will suffice 
to assume L violation only in the third generation, as a first approximation. 

The relevant contributions to these decays are given in Table 1. The Feynman dia- 
grams for the decays not involving taus, i.e. Xi ~^ ^sf f (/ = ^e, J^fi, u, d, c, s, b) 
are shown explicitly in Fig. 4.1. 



a) b) 




Figure 4.1: Feynman graphs for the decay Xi ^ ^3 f f where f ^ t. 

For this class of decays we have Z^, P°, and Sj exchange in the direct channel (Fig. 4.1a 
and b) and / exchange in the crossed channels (Fig. 4.1c and d). In particular in the 
case f = b the P° and Sj exchange contributions are significant. This is quite analogous 
to the results found in [185] for X2 ~^ Xiff decays. The particles exchanged in the s-, 
t-, and M-channel for the decays x? t^I'^i^u i} = Xi ~^ r^qq' {q,q' = u,d,s,c), 

Xi T^T^ui, and Xi ~^ Si^a are given in Tab. 4.1. 
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4.1 Decays of neutralinos 



Decay mode 


exchanged particle 


channel 




y qO pO 


s 




7 qO pO 


t 




z, sf, 


u 


X? I'l h {I = e, 


z 


s 




i>i 


t 






u 


Xi ^ / / (/ = e, u, d, s, c, b) 




s 




fl,2 


t 




/l,2 


u 


Xi ^3 


z, SI 


s 






t 






u 




w^, si 


s 






t 






u 


Xi^rqq' (q^u, c, q' = d, s) 


w^, si 


s 




Ql,2 


t 




91,2 


u 



Table 4.1: Contributions involved in the lightest neutrahno 3-body decay modes. The s-, 
t-, and ?i-channels are defined by: s — (pi — ^2)^, t — (pi — pa)^, and u — (pi — P4)^. See 
also Fig. 4.1. 
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In the calculations we have included all mixing effects, in particular the standard 
MSSM fi — fn mixing effects and those induced by the bilinear L violating terms, i.e. 
i?e(z>3) - /i° - /m(i>3) - ^0 - G^, [84], ff^^ - - [85], 1/3 - X° [186], and r - 
Xj mixings [58-60] . These mixing effects are particularly important in the calculations of 
the various L violating decay rates of x^, which are discussed in section 5.1.1. 



4.2 S quark decays 
4.2.1 Two body decays 

In Appendix B we give the Feynman rules for all vertices involving squarks, quarlcs, 
charginos and neutralinos. The decays of the six squarks are given both in the MSSM 
and the SSSM by 



g^X^/^{ml^,ml,mlo) 



167rm|^ 



A I. jki rjki 

- 4/iQ fn mqmpo 



m 



Ik 



m„ 



rripo 

3 



(4.1) 



167rm|^ 



+ ((/f)^ + (i7g) 



where Q — U,D refers to q and 

fjki 
JQ 

jjki 



-(v^Gg, + <^) 

Tjjki _ ^jki 



(4.2) 

(4.3) 
(4.4) 
(4.5) 



with G and H being either the SSSM or the MSSM couplings defined in Appendix B. 
4.2.2 Three body decays 

For definiteness and simplicity we assume only L violation in the third generation. A 

short discussion on ti ^ bl~^ in the three-generation model will be given at the end of 
Subsect. 5.2. In the one generation SSSM the stop mass matrix is given by 



^ {vuAt - novo + l^3V3) MIj + \vlht' + A[/R 



l.,2i 



with = \{g'^ - \g''^) {vl -'^1, + 'vf) and ^ur = \g'^{'^l -^1 + vj). The mass matrix 
for the sbottoms is given by 
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4.2 Squark decays 



1 ./2^ („,2 

^ 

eigenstates are obtained by {q = t,b): 



where Adl = -|(^^ + |^'^) (^^o - + ^oii = -]^f'^(^o - + ^D- The mass 
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The lighter stop three body channels are presented below. The complete formulas for 
the several contributions to each channel are given in [106] 
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(m6+m-o)2 



a 



167rm? sin'' 9w 



Jd s (Gx+x+ + <^x+t + G'x+b + + ^tb + Gil) 



(mi,+m-p)2 



a 



167rmf sin 9w 



Chapter 5 
Numerical Results 



5.1 Sugra Case 

For definiteness and simplicity we assume only L violation in the third generation case 
along this section. Note, in contrast, that in order to describe Flavor Changing Neutral 
Current (FCNC) effects such as the L conserving process — >^ cx? we need the three 
generations of quarks. 

The soft SUSY breaking parameters at the electroweak scale needed for the evalu- 
ation of the mass matrices and couplings are calculated by solving the renormalization 
group equations (RGB's) of the SSSM and imposing the radiative electroweak symmetry 
breaking condition. From the measured quark masses, CKM matrix elements and tan j3 
we first solve one-loop RGB's for the gauge and Yukawa couplings to calculate their cor- 
responding values at the unification scale. Assuming now universal soft supersymmetry 
breaking boundary conditions, we evolve downward the RGB's for all SSSM parameters, 
including full three-generation mixing in the RGB's for Yukawa coupling constants, as 
well as soft SUSY breaking parameters. Next, we evaluate the Higgs potential at the 
rut scale including the one-loop corrections induced by the Yukawa coupling constants of 
the third generation. The radiative electroweak symmetry breaking requirement fixes the 
magnitude of the SUSY Higgs mass parameter /xq and the soft SUSY breaking parameters 
Bq and B^. Notice that due to the minimization condition for in eq. (3.14) one can solve 
for as a function of ji-^. At this point, all L violating parameters at the electroweak scale 
are determined as functions of the input parameters (tan/9, mo, Aq, mi/2, sign(/xo), ""^t/s), 
where mo the common scalar mass and mi/2 is the common gaugino mass. Iteration is 
required because //o and are inputs to evaluate the loop-corrected minimum. Having 
determined all parameters at the electroweak scale, we obtain the masses and the mixings 
of all the SUSY particles by diagonalizing the corresponding mass matrices. At this stage 
we also choose /is in order to get a sufficiently light v^. 

5.1.1 Neutralinos 

In this section we present numerical predictions for the lightest and second lightest neu- 
tralino production cross sections in e+e~ collisions, namely, e+e~ XiXi, XiX2- Moreover 
we will characterize in detail all branching ratios for the hghtest neutralino decays in the 
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5.1 Sugra Case 



SSSM. It will suffice to assume L violation only in the third generation, as a first approx- 
imation in this case. 

The relevant parameters include the L violating parameters and the standard Minimal 
SUGRA (MSUGRA) parameters mi/2, rno, tan/?. The absolute value of /iq is fixed by 
radiative breaking of electroweak symmetry. We take positive to be in agreement with 
the b ^ s'j decay [91]. As representative values of tan/3 we take tan/? = 3 and 50. It is a 
feature of models with purely spontaneous breaking of R-parity that neutrinos acquire a 
mass only due to the violation of R-parity [17, 22, 186]. This feature also applies to the 
SSSM. As a result the L violating parameters are directly related with m,^^, the mass of 
the neutrino u^, which is generated due to the mixing implicit in (3.4). 

Neutralino Production 

While the L violation would allow for the single production of supcrsymmctric parti- 
cles [58-60], for the assumed values of the L violation parameters indicated by the sim- 
plest interpretation of solar and atmospheric neutrino data [26, 144-146], these cross 
sections are typically too small to be observable. As a result neutralino production at 
LEP2 in our model typically occurs in pairs with essentially the same cross sections as in 
the MSUGRA case. In Fig. 5.1a and b we show the maximum and minimum attainable 
values for the e~^e~ XiXi ^^'^ — * X1X2 production cross sections as a function 
of m^o at = 205 GcV. We compare the cases tan/3 = 3 and tan/5 = 50, varying 
mi/2 between 90 GeV and 260 GeV and mo between 50 GeV and 500 GeV. One can see 
that, indeed, these results are identical to those obtained in the MSUGRA. The XiXi 
production cross section can reach approximately 1 pb. In our calculation we have used 
the formula as given in [185] and, in addition, we have included initial state radiation 
(ISR) using the formula given in [187]. Note that and are exchanged in the t- and 
M-channel implying that a large fraction of the neutralinos will be produced in the forward 
and backward directions. 

In order to show more explicitly the dependence of the cross sections on the parameters 
mo and mi/2 we plot in Fig. 5.2a and b the contour lines of a{e^e~ XiXi) iii the mo- 
mi/2 plane at a/s = 205 GeV for tan/5 = 3 and tan/3 = 50. The contour lines for 
(j(e"'"e~ — >• X1X2) given in Fig. 5.2c and d. 

Neutralino Decay Length 

For the observabihty of the L violating effects it is crucial that with this choice of pa- 
rameters the LSP will decay most of the time inside the detector. The neutrahno decay 
path expected at LEP2 depends crucially on the values of L violating parameters or, 
equivalently, on the value of the heaviest neutrino mass, m,^.^. We fix the value of m^^g as 
indicated by the analysis of the atmospheric neutrino data [26, 144]. It is important to 
note that, as explained in [34], due to the projective nature of the neutrino mass matrix 
[22], only one of the three neutrinos picks up a mass in tree approximation. This means 
that, neglecting radiative corrections which give small masses to the first two neutrinos 
in order to account for the solar neutrino data, the neutralino decay length scale is set 
mainly by the tree-level value of m^^ . In ref . [34] we have explicitly shown that this is a 
good approximation for most points in parameter space. 
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a) a [fb] b) ^ ^£5] 




40 50 60 70 80 90 100 40 50 60 70 80 90 100 

m^o [GeV] m^^o [GeV] 

Figure 5.1: Maximum and minimum attainable values for the e+e" XiXi (f^H lines) 
and e'^e~ X1X2 (dashed lines) production cross sections in fb as a function of m^o for 
= 205 GeV, 50 GeV < mo < 500 GeV, 90 GeV < mi/2 < 270 GeV, a) tan (5 = 3,'and 
b) tan /3 = 50. ISR corrections are included. 
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Figure 5.2: Contour lines of the production cross sections in fb, in the mo-mi/2 plane for 
= 205 GeV, a) e+e" ^ xSx?, tan/? = 3, b) e+e" ^ xSx?, tan/5 = 50, c) e+e" ^ xlxl 
tan/3 = 3, and d) e'^e" XiX2: tan/3 = 50. ISR corrections are included. 
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a) L(x?) [cm] b) ^^(^O) [^^j 




Figure 5.3: Decay length of the hghtcst neutrahno in cm for = 205 GeV, a) as a 
function of rrii,^ for m^o = 60, 70, 80, and 90 GeV, b) as a function of m^o for 777.1,3 — 
0.01,0.1, and 1 eV. 



In Fig. 5.3 we plot the Xi decay length in cm expected at LEP2 for ^/s — 205 GeV. 
Here and later on we consider the neutralinos stemming from the process e'^e" — > xlxl 
when discussing the decay length. In Fig. 5.3a wc plot the Xi decay length in cm as a 
function of neutrino mass mi^g, for different 777^0 between 60 and 90 GeV, with ?77o = 100 
GeV, and tan/5 = 3. As can be seen the expected neutralino decay length is typically 
such that the decays occur inside the detector, leading to a drastic modification of the 
MSUGRA signals. An equivalent way of presenting the neutralino decay path at LEP2 
is displayed in Fig. 5.3b, which gives the decay length of Xi as a function of ?77^o for 
777,^3 = 0.01, 0.1, and 1 eV. Finally, we show the dependence of the neutralino decay 
path on the supergravity parameters fixing the magnitude of L violating parameters or, 
equivalently, the magnitude of the heaviest neutrino mass, 7772,3. ^^S- 5.4a and b we 
plot the contour lines of the decay length of x1 in the ■mo-'mi/2 plane for m^^ — 0.06 eV, 
tan /3 = 3 and 50. Note that the decay length is short enough that it may happen inside 
typical high energy collider detectors even for the small neutrino mass values ~ 0.06 eV 
indicated by the atmospheric neutrino data [26, 144]. For large values of tan/3 the total 
decay width increases and, correspondingly, the decay path decreases due to the tau 
Yukawa coupling and the bottom Yukawa coupling. 

Neutralino Branching Ratios 

As discussed in the beginning of this section, the lightest neutralino Xi will typically 
decay in the detector. In the following we present our results for the branching ratios of 
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Figure 5.4: Decay length of the hghtest neutrahno in cm in the mo-mi/2 plane for sfs = 
205 GeV, a) tan/3 = 3, and b) tan/? — 50. The L violating parameters are fixed such 
that = 0.06 GeV. 



all L violating 3-body decay of Xi, and of the radiative decay Xi ~^ ^37- The Feynman 
diagrams for the decays x5 ~^ ^sff if — 6, J^e, l-i', u, d, c, s, b) are shown in Fig. 4.1. 

In the following plots Fig. 5.5 - 5.12 we show contour lines in the mo-mi/2 plane for 
the branching ratios in % of the various Xi decays, in (a) for tan/? = 3 and in (b) for 
tan/? = 50. We have fixed the mass of the heaviest neutrino to rrii,^ = 0.06 eV [26, 144]. 
It turns out, that in the range 10~^ eV < m^^g < 1 keV all the Xi decay branching ratios 
are rather insensitive to the actual value of m^^. This is an important feature of our 
supergravity-typc SSSM. It is a consequence of the fact that, as a result of the universal 
supergravity boundary conditions on the soft breaking terms, all L violating couplings are 
proportional to a unique common parameter which may be taken as Hs/no- For a more 
detailed discussion on this proportionahty the reader is referred to ref. [34]. Also note 
that for mi/2 ^ 220 GeV the neutralino mass becomes larger than mw and niz so that 
Xi decays into real W and Z are possible. The effects of these real decays can be seen 
for mi/2 ^ 220 GeV in most of the following plots. For the large tan/? case (tan/? = 50) 
and mi/2 ^ mo the mass of the lighter charged boson is smaller than m^o (upper left 
corner of Fig. 5.5b - 5.12b). In this region of the parameter space the two 2-body decays 
Xi — > W^T^ and Xi ~^ ^t"^^ compete. The first one is L violating, but has more phase 
space than the second one which is L conserving, since is mainly a stau. For this 
reason, the most import final state is t~^t~i'3, followed by r^qq' and r^l^Ui {I = e,fi) as 
shown in Figs. 5.11, 5.9, and 5.8, respectively. All other final states have nearly vanishing 
branching ratios in this corner of the parameter space. 
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a) mi/2 [GeV] b) m^i^ [GeV] 




100 200 300 400 500 100 200 300 400 500 

mo [GeV] mo [GeV] 

Figure 5.5: Branching ratios for x5 in % in the mo-mi/2 plane for a) tan/3 = 3, 

and b) tan /3 = 50. The L violating parameters are fixed such that vriv^ — 0.06 GeV. 

Fig. 5.5a and b exhibit the contour lines for the branching ratio of the invisible decay 
X? — 3 z/. This branching ratio can reach 7% for the parameters chosen. In Figs. 5.6, 

5.7 we show the branching ratio for the decays Xi ~^ ^z^^ ^~ and Xi ^ ^z^l^l where / 
and q denote the leptons and quarks of the first two generations, summed over all flavors. 
These branching ratios can go up to 3% and 15%, respectively. Notice that the sneutrino, 
slepton, and squark exchange contributions to the x5 decays become larger with increasing 
mo, despite the fact that the increase of the scalar masses m^, m^-, rriq suppresses these 
exchange contributions. This trend can also be observed in Fig. 5.7, 5.8 and 5.9. This 
happens because the tadpole equations correlate /iq to mg. Increasing /xq while keeping 
Ml and M2 fixed implies increasing the gaugino content of and, hence, enhancing the 
X?-/-/ couphngs. 

In Figs. 5.8, 5.9 we show the contour lines for the branching ratios of the LSP decays 
involving a single tau, namely Xi ~^ ^it^P and Xi 't^QQ', where I, q, and q' are 
summed over the first two generations. The branching for these decay modes can reach 
up to 20% and 60% respectively. For mi/2 ^ 220 GeV decays into real dominate. If 
this is the case and if both x? produced in e+e~ XiXi decay according to these modes 
this would lead to very distinctive final states, such as Ajt'^t'^, t'^t'^1~1~ {I — e,/i), 
or r~^T~^e~ n~ . The full list of expected signals is given in Tab. 5.1. The first column 
in this table specifies the two pairs of x? decay modes, while the second one gives the 
corresponding signature. In the last column we state whether the corresponding signature 
exists for e+e" — > X1X2 production within MSUGRA. 
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a) mi/2 [GeV] b) ^1/2 [GeV] 
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mo [GeV] mo [GeV] 

Figure 5.6: Branching ratios for Xi ~^ ^~ in % in the mo-mi/2 plane for a) tan /9 = 3, 
and b) tan P = 50. Here / is the sum of e and /i. The L violating parameters are fixed 
such that m,y, = 0.06 GeV. 
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a) mi/2 [GeV] b) m^i^ [GeV] 




100 200 300 400 500 100 200 300 400 500 



mo [GeV] mo [GeV] 

Figure 5.7: Branching ratios for y^^^ v-iqq\Ti.%\VL the mo-mi/2 plane for a) tan/3 = 3, 
and b) tan /3 = 50. Here q is the sum over d, s, and c. The L violating parameters are 
fixed such that m,^, = 0.06 GeV. 



The LSP decays involving only third generation fermions, namely, Xi ^ ^3 b b and 
Xi 1^3 T~ are different from those into the first and second generation fcrmion pairs, 
because the Higgs boson exchanges and the Yukawa terms play a very important role. 
This can be seen in Figs. 5.10, 5.11 , where we plot the contour lines for these decays. The 
branching ratio of x5 ^^3 can reach up to 97%. The decay rate is large because the 
scalar exchange contributions (5"°, Pj^, hk) are large for mi/2 200 GeV. Note that this is 
also the case for tan /9 = 3, because not only the neutrino-ncutralino mixing proportional 
to iTLy.:^ is important but also the neutrino-higgsino mixing proportional to fis/fiQ. The 
decrease of the branching ratio with increasing mo is due to the decrease of the higgsino 
component of Xi the increase of the Higgs boson masses. For mi/2 ^ 200 GeV the 
decays into real and Z° are possible, reducing the branching ratio of X? — > t's b. 
As shown in Fig. 5.11 the branching ratio for Xi ~^ ^3 t~ is very small for tan/? = 3 
and mi/2 ^ 200 GeV. This is due to the destructive interference between contribution 
and the contributions of the exchanged charged scalar particles (mainly due to the stau 
components of S^). 

Finally we have also considered the radiative LSP decay mode x5 ~^ ^3 7 [93]. In 
Fig. 5.12 the branching ratio for this mode is shown. This decay proceeds only at one- 
loop level and therefore is in general suppressed compared to the three-body decay modes. 
However, for mi/2 % 125 GcV and large tan/3 it exceeds 1%, leading to interesting 
signatures like e+e~ — > XiXi ~^ f^^ ^ + /r- Due to initial state radiation it can 
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a) mi/2 [GeV] b) ^1/2 [GeV] 
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mo [GeV] mo [GeV] 

Figure 5.8: Branching ratios for Xi ~^ J^i P in % in the mo-mi/2 plane for a) tan /9 = 3, 
and b) tan (3 = 50. Here / is the sum of e and /i. The L violating parameters are fixed 
such that = 0.06 GeV. 
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a) mi/2 [GeV] b) ^1/2 [GeV] 




100 200 300 400 500 100 200 300 400 500 



mo [GeV] mo [GeV] 

Figure 5.9: Branching ratios for — > g^' in % in the mo-mi/2 plane for a) tan/3 = 3, 
and b) tan P = 50. Here q is the sum over u, d, s, and c. The L violating parameters are 
fixed such that m^, = 0.06 GeV. 
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a) mi/2 [GeV] b) [GeV] 




100 200 300 400 500 100 200 300 400 500 

mo [GeV] mo [GeV] 

Figure 5.10: Branching ratios for ^ 1/3 bb in % in the too~"^i/2 plane for a) tan (3 — 3, 
and b) tan /3 = 50. The L violating parameters are fixed such that m,^^ — 0.06 GeV. 

easily happen that a second photon is observed in the same event. 

Signals at e"^ e~ colliders 

The complete list of possible signatures stemming from LSP decays in SSSM is shown 
in Tab. 5.1. In this table we also indicate whether the same signatures could also arise 
in MSUGRA as a result of e+e~ X1X2 followed by the MSSM decay modes of if 
its production is kinematically allowed. The final states 4 jets + [188], r + 2 jets + 
^T, and r + (e or /i) + would also occur in MSUGRA via the decay of X2 ^^^'^ xf- 
However one expects in general that these decay modes are suppressed within MSUGRA. 
In contrast in the SSSM these signatures can be rather large as can be seen from Figs. 5.8, 
5.9. Note moreover, that some of the L signatures are practically background free. For 
example, due to the Majorana nature of x1, one can have two same-sign r leptons + 4 
jets + ^T- Other interesting signals are: r + 3 (e and/or /i) + ^t, 3 r + (e or /x) + ^t, 
T + (e or /x) + 2 jets + fr, r + A jets + ^t, t^t^ + (e or /x) + 2 jets + ^t, or t^t^ 
+ r^l'^ + with I — e,iJ,. In Table 5.2 we give masses and branching ratios for typical 
examples. 

As it is well known, also in gauge mediated supersymmetry breaking models (GMSB) 
[66] the neutralino can decay inside the detector, because the gravitino G is the LSP. It 
is therefore an interesting question if the SSSM can be confused with GMSB. To answer 
this question let us have a look at the dominant decay modes of the hghtest neutrahno in 
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Combination or xi decay modes 


signature 


MSUGRA-like 


(3i/) (3 1/) 




yes 


(3i/) {U3I+I ) 


2 leptons + 


yes 


(3i/) (i/3?g) 


2 jets + 


yes 








(3 i/j [I'lT with I — e, ijl 


T + (e or //) + 


yes, but suppressed 


(3 v) {r^qq') 


T + 2 jets + 


yes, but suppressed 


(3z/) (z/3 7) 




yes 


/ 7-1-7 \ / 7/-I-7/ — \ 

(2/3/+/ ) {p^l'^V ) 


4 leptons + 


no 


( 7-1-7 \ / \ 

(2/3/+/ ) [y-iqq) 


2 leptons + 2jets + 


no 


( 7_|_ 7 \ / 7 T\ 






yV-^l I ) [yiT ) Wltil L — e, /i 


T + 3 (e and/or /i) + 


no 


{jy3T~^T~) (uiT^P) with I = e, n 


d r + (e or /x) + 


no 


{r^qq') 


r + 2 leptons + 2 jets + 


no 


(z/3r+r ) (r^gg ) 


3 r + 2 jets + IfT 


no 


(7/3/+/-) (7/37) 


2 leptons + 7 + -^t 


no 


/ — \ / — \ 








4 jets + 


yes, but suppressed 


(1/300) (t'sOOj 






[usqq) [uiT^l+) with I = e, /i 


T + (e or //) + 2 jets + 


no 


(1/300) \yiT l^) with I — e, ijl 








T + 4 jets + 


no 


(1/36&) (T=^gg) 








2 jets + 7 + 


no 


(1/3&6) (1/37) 








_it_± 1 1^11'^ 1 

T T + /^/ + 


no 






no 


{vitH^) {r^qq') 


T=^T=^ + (e or //) + 2 jets + 


no 




T^T^ + (e or //) + 2 jets + |V 


no 




T + (e or //) + 7 + 


no 


{r^qq') {r^qq') 


T=^T=^ + 4 jets + 


no 




t^t"^ + 4 jets + 


no 


{r^qq') (z/37) 


r + 2 jets + 7 + 


no 


(^^37) (^^37) 


27 + ^T 


no 



Table 5.1: The signatures expected from the process e+e — > XiXi iii the SSSM. 
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Figure 5.11: Branching ratios for Xi — > ^3 t+ t~ in % in the mo-mi/2 plane for a) tan (3 ■ 
3, and b) tan (3 — 50. The L violating parameters are fixed such that m,^^ — 0.06 GeV. 



GMSB. If the hghtest neutrahno is the NLSP, its main decay mode in GMSB is 



where G is the gravitino. For the case where at least one of the sleptons is lighter than 
the hghtest neutralino the latter has the following decay chain x\ l"^ l^V^G 
In principle three-body decay modes mediated by virtual photon, virtual Z-boson and 
virtual sfermions also exist. However, in the neutralino mass range considered here these 
decays are phase-space-suppressed [66, 189]. This implies that the SSSM can not be 
confused with GMSB, because (i) in GMSB the final states containing quarks are strongly 
suppressed, and (ii) GMSB have lepton flavor conservation, and therefore there are no final 
states like e+e+r^r" + ^t- A further interesting question would be how the neutralino 
phenomenology changes in a GMSB scenario with L violation. The main consequence 
would be an enhancement of final states containing photons and/or leptons. A detailed 
study of this question is, however, beyond the scope of the present work. 

5.1.2 Lightest Stop Two-Body Decays in SUGRA 

For definiteness and simplicity we assume only L violation in the third generation case. 
The discussion on t\ — > hl^ in the three generation model will be given at the end of 
Subsect. 5.2 

Note, in contrast, that in order to describe Flavor Changing Neutral Current (FCNC) 
effects such as the L conserving process ti ^ c x° we need the three generations of quarks. 
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Figure 5.12: Branching ratios for x° — > 1/37 in % in the mo-mi/2 plane for a) tan/3 = 3, 
and b) tan /3 = 50. The L violating parameters are fixed such that — 0.06 GeV. 



We scan the soft SUSY breaking parameter space in the range 

mo < 700 GeV 
100 GeV < mi/2 < 400 GeV 

\Aq\ < 1000 GeV, (5.1) 
m^3 < 18 MeV 
1.8 < tan/3 < 60 

the previous range on tan /? guarantee that both hf and hy, will be perturbative. For the 
CKM matrix, we use the Particle Data Group convention [190], taking K^s = 0.2205, 
Kcb — 0.041, \Kub/Kcb\ — 0.08 and neglecting CP violation, i.e. 5 — 0. Notice that here 
we scan over a much larger range for /i^ than used in ref. [70]. 

The resulting region of lightest stop and chargino masses is displayed in Fig. 5.13. 
Neglecting the three-body decays, we find that in Region 1 of the mf^-m-+ plane, BR{ti — > 
cx?) + BR{h br) ^ 1. In Region II BR{ii ^ br) + BR{ii bxf) « 1 (i=l,2). In 
Region III BR{ii bT)+BR{ti b + xj)+BR{ti ^tv^)^l (i=l,2), while in region 
IV BR{h ^bT)+ BR{ii ^ bxt) + BR{h ^tiys) + BR{h ^ i x°) ~ 1 (j = 1, ■ ■ ■ , 4). 
Note that in each region the exact equality to 1 is reached when the FCNC processes are 
fully included. 

In Sec. 4.2 we give the two-body squark decay-widths, for squarks of all three genera- 
tions. These equations reduce to the expressions found in Ref. [77] provided one identifies 
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tan /3 = 


q 


tan /3 = 50 




A 


B 


\j 


A 


B 


c 




54.6 


59.0 


no 


60.0 


61.5 


94.4 


777.50 


91.0 


96.8 


102.9 


107.2 


111.1 


116.4 




180.5 


449.6 


466.0 


178.2 


448.7 


465.1 


men 


170.5 


445.7 


4o(J.D 


171.6 


446.1 


451.0 


rrie^ 


194.2 


455.2 


471.4 


195.3 


455.8 


471.9 


ruq 


398.1 


572.8 


TACT y< 

705.4 


398.1 


572.8 


705.4 




261.4 


328.5 


44z.z 


279.9 


355.2 


466.3 




361.3 


479.1 


61z.l 


243.0 


343.0 


470.1 


7-) 7->/ ~n \ 


0.5 


4.5 


1.8 


0.3 


1.2 


1.9 




0.2 


1.1 


0.5 


0.2 


0.3 


0.6 




1.0 


8.6 


4.0 


0.5 


2.2 


4.4 




0.6 


5.6 


18.0 


0.5 


1.8 


17.8 




1.1 


16.1 


53.7 


0.9 


5.1 


53.2 




96.5 


62.6 


13.4 


97.1 


88.4 


13.3 




0.1 


1.5 


8.6 


0.5 


1.0 


8.8 



Table 5.2: Masses and branching ratios for the points: A (mi/2,'^o) = (153,155), B 
(mi/2,mo) = (153,440), and C (mi/2,mo) = (251,440) for both tan/3 = 3 and 50. The 
masses are given in GeV and the branching ratios in % and we only give those larger than 
0.1%. Here the same summations of the final states are performed as in the figures, rriq 
is the averaged squark mass for the first two generations. 



= "^o^ ^iid r^^^ = sin Oi. They also generalize the results for the SSSM to the 
three-generation case. 

In the MSSM the main ii decay channel expected in region I of Fig. 5.13 is the loop- 
induced and flavor-changing — > c [71-75] . As is well-known, the FCNC processes in 
the MSSM in general involve a very large number of input parameters. For this reason, 
following common practice, we prefer to perform the phenomenological study of flavor 
changing processes in the framework of a supergravity theory with universal supersym- 
metry breaking. The simplest description of FCNC processes in SUGRA models uses the 
so-called one-step approximation. Here we start by reproducing the standard calculation 
for cx\ as in [71]. To do this consider only the effect of the third generation Yukawa 
coupling. From our general eq. (4.1) we have for ti = ui 



,2 



nh^cx'D^^iTuLi.f 



mr 1 - 



w 



(5.2) 
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Figure 5.13: Kinematical regions in the mr -m-+ plane. From left to right: Region I mz < 
m-+ + rrif,; Region II m-+ + m;, < mr < mt, Region III mt < mr < m^o + mu and region IV 
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m^^ > m^o + mt 



with 



UL13 



Al cos 9^ — Aji sin 

2 2 

— mj 



In the one-step approximation A/,, A/j are given by 

tu 



Al = (M, 



C/;23 ~ (M2)23 



U)2& ~ (^L/)23 



167r2 



15^2— - ---- ■ 2- 



(5.3) 

(5.4) 
(5.5) 



with tu = ln{MG/mt). So, in the one-step approximation we have 



r(ti ^ c x^) ^ Fhti6^2 COS Oi - 5a sin Oi) 



V2 



(tan 61^^ A^ii + 3iV; 



12 



Xi_ 



(5.6) 



where the pre-f actor F = (^^^KcbKtb) ~ 6 x 10 and the parameter 5^2 is given 
by 



Om2 = K K ~ i 



2 2 

mi — mf 



(5.7) 



is basically independent of the initial conditions due to the mo dependence both in the 
numerator as in the denominator and 



ml — m^ 

cl 



tL 
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Note however, that the one-step approximation includes only the third generation Yukawa 
couplings and neglects the running of the soft breaking terms [71-76]. Such an approx- 
imation is rather poor for our purposes, since we will be interested in comparing with 
L violating decay modes (see section 5.1.3). In order to have an accurate calculation of 
the respective branching ratios we need to go beyond the one-step approximation. We 
therefore use a exact numerical calculation for the FCNC process t — ex? in which the 
running of the Yukawa couplings and soft breaking terms is taken into account. First 
we have checked that indeed the effect of the Yukawas from the two first generations is 
negligible. However the same is not true for the running of the soft breaking terms. As 
can be seen from Figure 5.14 the range of variation that we obtain from the numerical 
solution is 

r(fi^cx?)~(10-'^-10-^)GeV (5.9) 

depending on the assumed value of m\/2 and tan (5. In this figure we have compared the 
decay width obtained from eq. (4.1) with the approximate formula in eq. (5.2) for two fixed 
values of mi/2, tan/? and taking Aq = 0. The approximate formula only reproduce well 
the numerical result for the academic case of no SUSY breaking gaugino mass, mi/2 = 0. 
For the reahstic case mi/2 > 100 GeV, the exact solution is usually one decade smaller 
than the approximate one. In the one-step approximation r{ii — > cXi) can be arbitrarily 
small if the two terms 5^2 cos % and S^sinOi in eq. (5.6) cancel. This behavior can be 
illustrated in Figure 5.14 by the dashed line labeled 358, which corresponds to mo = 358 
GeV. One sees clearly that while the approximate solution goes to zero, the numerical 
one reaches a minimum value around 10~^^ GeV. The wrong behavior of the approximate 
solution indicates that the Sa depends strongly on the scale. For example, the RGE for Ai, 
is very sensitive on mi/2 and tan/3 and in the one-step approximation there is no explicit 
dependence on mi/2, which is crucial. Both solutions increase with tan (3, as expected by 
the bottom Yukawa dependence explicit in eq. (5.6) and remain practically constant for 
large mi/2 values. 

5.1.3 Two-Body Decays of the Lightest Stop: the SSSM case 

In contrast to the case of an L conserving supergravity theory (MSUGRA), in the L 
violating through Superrenormalizable terms (SSUGRA) case, one can have a competing 
L violating stop decay mode in region I of Fig. 5.13. From eq. (4.2) with r = F3" one can 
easily compute the L violating stop decay width ii — > br, 

r{ti^br) = I6nm^ {-4=U^2hbC0.{V^2htS0- - V^^C0-)mbmr 

h 

+ m2htSo, - V,\cof + U;fhlc%]iml ~ ml ~ ml)} (5.10) 

which coincides with the result found in Ref. [70]. In [107] it was shown that, except for 
Us2 which determines the SU(2)-conserving mixing of the Higgsino with the left-handed r, 
all other mixing matrix elements V^i and U^i are proportional to v^^ (defined in eq. (1.14) 
)and therefore to the mass. Neglecting these terms we have from eq. (5.10) 



777^ TTl^ TTl^ j 

r(ti - 6r) ^ - sin^ i'hlcUml -ml- ml) (5.11) 



tl 
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Figure 5.14: Comparison between the exact numerical calculation (ordinate) and the one-step 
approximation (abscissa) for the ii — > ex? decay width for various values of tan/3 and mi/2 
with Aq = and niQ varying in the indicated range. The dotted left diagonal line would signify 
the equality between the estimates, while the right diagonal line would indicate one order of 
magnitude difference. Results of both estimates indicated in the lower right legend. More 
details are found in the text. 



noting that, to a good approximation, 

sin^'l (5.12) 

where /is corresponds to the superrenormahzable parameter in basis I defined in sec- 
tion 1.2. The lesson here is that the L violating decay rate T{ti -^hr) is proportional to 
/i3 or, equivalently, to sin^^', instead of sin^^, and thus not necessarily small, since it is 
not directly controlled by the neutrino mass. In other words, there can be cancellations in 
the latter but not in the L violating branching ratio. 

The meaning of the factor sin^ ^'hl may also be seen in basis II, where e^^ = 0. In this 
case is proportional to the mass so that, as already mentioned, in this basis all the 
elements U^i and Vsi are small [92]. Neglecting these terms, T{ti br) may be written 
directly from the interaction term tLbRTi, which is induced by the trilinear term in the 

= 0-basis given in eq. (1.16) as 

^¥ = h = hsin^' (5.13) 

which is the factor in eq. (5.11). Note, however, that in our numerical calculation to be 
described below we have used for r(ti br) the full expression given in eq. (4.1). 

Next section we will determine the conditions under which the L violating decay width 
Till br) can be dominant over the L conserving ones, T(ti cXi) and T(ti bxf). 
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(5.14) 



Using the eq. (5.11) and neglecting charm, tau and bottom masses we get 




2 



- 10 



,-5 



sin^ ^' 



(5.15) 



Therefore T(ii ^ c x^) wiU start to compete with T{ii br) from sin ^' ^ 5 x 10~^ (10~^) 
for tan/5 large (small). In Fig. 5.15 we compare BR{ti — > cXi) [calculated numerically 
from their exact formula in (4.1)] with BR{ii br) within the restricted region of the 
mfj-m^o plane where only those two decay modes are open. We consider different 171^^ 
values (these correspond to relatively small values of the L violating parameters lusl, {v^l < 
1 GcV). We vary the SSSM parameters randomly obeying the condition m^^ < m-± 
and depict the corresponding region in light grey. The upper-left triangular region is 
defined by kinematics and corresponds to m^^ < m^o +mc, so that BR{ti br) — 100%. 
The lower-right grey corresponds to m^^ > m^o + rric when the sampling is done over 
the region defined by eq. (5.1). One notices from Fig. 5.15 that in the central region 
the dominant stop decay mode is ii br with branching ratio BR(ti ^ br) > 0.9. 
The dotted lines in the light grey region indicate maximum 1/3 mass values obtained in 
the scan. In the calculation of the 1/3 mass, we have allowed only up to one order of 
magnitude of cancellation between the two terms which contribute to sin^. Therefore if 
the lightest stop only decays into the two modes considered here, the processes fi — > &t, 
will be important even for the case of very light z/3 masses. 

We note however that we can use the limits obtained from leptoquark searches [191] 
in order to derive limits on the top-squark in the SSSM. In particular, if BR{ii — > br) — 1 
stop masses less than 99 GeV are excluded at 95% of CL., under the assumption that 
the three-body decays of the stops are negligible. Therefore, the dark region in Fig. 5.15 
would be ruled out. In ref. [192] we have determined the corresponding restrictions on 
the SUGRA parameter space. 

The dependence on the 1/3 mass may be seen in Fig. 5.16 where the role played by tan j3 
is manifest. In this figure we have shown BR{ii — > 6t) as function of the lighter stop 
mass for z/3 mass in the sub-cV range, indicated by the simplest oscillation interpretation 
of the Supcr-Kamiokande atmospheric neutrino data. We have obtained such z/3 mass 
values numerically, allowing only one decade of cancellation between the two terms that 
contribute to sin^ in eq. (1.26). The degree of suppression for Hs/no obtained numerically 
agrees very well with the expectations from the approximate formula for the minimal 1/3 
mass in eq. (1.29). In contrast with Ref. [70], in our case BR{ti — > br) decreases with 
tan/3. The reason for this difference is that here wc take into account the fact that the 
mixing parameter Tuin obtained from the RGE depends on hi in eq. (5.2), while in 
ref. [70] was simply regarded as a phenomenological input parameter (called 5 there). 
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Figure 5.15: Regions where the l\ ^ hr decay 
branching ratio exceeds 90% in the mi^-rfiyQ^ plane 
for different niv^ values. The SSSM parameters 
are randomly varied as indicated in the text under 
the restriction mr < to-± + to?,. The upper-left 
triangular region corresponds to m^^ < m^o + ra^ 
so that only the ti — > &t decay channel is open. 
The lower-right unshaded region corresponds to 
mr > TO~+ + mf,. 



Figure 5.16: J5i?(ii br) as function of the 
lighter stop mass for 1^3 mass in the sub-eV range 
and two different values of tan/3 and /LtsZ/io- This 
prediction is natural in the sense that we have al- 
lowed only up to one order of magnitude of can- 
cellation between the two terms that contribute to 
sin^. 



The message from this subsection is that in our SSUGRA the L violating decay mode 

.15 



ti — ^ 6r can very easily dominate the L conserving decay mode ti —>■ cXi, even for very 



small neutrino masses. 
Region II 

In region II the L conserving decay mode ii —>■ bxt is open (but not ti —>■ tu), and 
competes with the L violating mode ii br. Replacing the subindex 3 by 1 on the 
diagonalization matrices U and V in eq. (5.10) we get the corresponding expression for 
T{ti —* bxt)- In order to get an approximate expression for the ratio of the two main 
decay rates in this region, we note that in SUGRA with universality at the unification 
scale, the lightest chargino is usually gaugino-like, implying that ~ 1. In addition, 
the lightest stop is usually right-handed, hence sin^ 9^ > cos^ 9^. This way we find 

r(r) _ Tih^br) ^ sm'^'hlcos'9i ^ ^^^^^^ 



T{xt) m^bxt 



{V{-^ cos 9i - V^Jit sin 9iy + Ulfhl cos^ 9i 



where is a kinematical factor depending on the lightest stop and chargino masses, and 
here we have defined ht^b = ht,b/ g- The presence of the bottom quark Yukawa coupling 
indicates that large values of tan (3 are necessary to have large L violating branching ratios 
in this region. In fact, we have checked numerically with the exact expressions that in 
Region II (RII) r(r)/r(xi') ^ 1 only for large tan/3 as we will see in the next figures. 
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In Fig. 5.17 we show the regions in the m^±—mi^ plane wliere BR{ti ^ br) dominates 

over BR{ti bxt)- In the upper-left region the decay mode ii — > bxt is not allowed and 
corresponds to Region I. Below and to the right of this zone, and above and to the left 
of three rising lines, lies region RII where r(r)/r(x^) > 1. The three lines correspond 
to I /is I < 80 GeV (dashed), I/is] < 60 GeV (dotted), and I/is] < 40 GeV (dot-dashed), 
respectively. The proximity to the upper-left zone indicates that the L violating decay 
dominates only close to the threshold where there is a high kinematical suppression of 
the L conserving one, through the factor K. Unlike the case of region I this requires 
large values of the L violating parameters. Note, moreover, that if the stops have a small 
mixing {cos9t ^ 0), then T{t)/T{xi) < 1 in RII. 



TRTp 




W<40GeV 



eo 90 100 110 120 130 140 ISO 160 170 IflO 
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Figure 5.17: Contours of BR{ti br) > BR{ii bxt) in the nir -m-+ plane for {v^l < 10 
GeV. Three different maximum values for {fi^l are considered: l/^al < 40 GeV (dot-dash), {fi^l < 
60 GeV (dots), and I/U3I < 80 GeV (dashes). The region where mi < rrih + m~+ corresponds to 
the previously studied Region I. 

A simpler expression for the ratio of decay rates in eq. (5.16) is obtained if we take 
Vii ~ 1 and assume no kinematical suppression in eq. (5.16) through the factor K: 

^^^^ sin'fhl (5.17) 



Note that the presence of the parameter sin^' = /iaZ/i indicates that the L violating 
decay mode is not strictly proportional to the neutrino mass, but proportional to the 
SSSM parameter 

However generically we expect some correlation with the z/3 mass, especially in the 
case where the boundary conditions in the RGE are universal and there are no strong 
cancellations between two terms that contribute to sin^ as shown in Fig. 5.18. In this 
figure we plot the ratio T{t) /T{xi) in RII as a function of the z/3 mass. Both decay rates 
have been calculated numerically from the exact formulas. In this figure we have imposed 
both = M| and Bq = B^ within 0.1% at the GUT scale. Cancellation between 
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the Am^ and Ai? terms in the neutrino mass formula of eq. (1.26) are accepted only 
within 1 decade. As a reference we have drawn the line corresponding to AB = and 
Am^ = Am^j„ (Am^ is negative and its magnitude is bounded from below by Am^j„) 
at the weak scale, which gives an idea of the value of the neutrino mass when there is no 
cancellation between the AB and Am^ terms. 

We have imposed an upper bound on m^,^ at the collider experimental limit of the 1/3 
mass, and have chosen fixed values of /isZ/io = 1, 0.1, and 0.01. The allowed region for 
/^s/z^o = 1 is above the dashed line. In the case of fis/fio = 0.1 (0.01) the allowed region 
lies enclosed between the solid (dotted) lines. The effect of tan P is to increase the ratio 
r(r)/r(x^): the minimum value of the ratio is obtained for tan/? ^ 2 and the maximum 
corresponds to tan f3 ^ 60. The extreme values of tan/5 are dictated by perturbativity. 
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Figure 5.18: Regions for r(ti hT)/T{ii 
b xt) ^•s a function of the ^3 mass with the uni- 
versality condition Bq = B3 at the unification 
scale imposed at the 0.1% level as indicated. 
Its effect is to alter the maximum attainable 
f3 mass. The dot-dashed line corresponds to 
the case where AB = at the weak scale. 



Figure 5.19: Regions for r{ii 6r)/r(ti 
bxi) as a function of the u-^ mass for different 
levels of cancellation between the two terms 
that contribute to the neutrino mass. We im- 
pose the universality condition m^^ = M£ at 
the unification scale, but Bq and B^ are not 
universal. We take |/i3//io| = 1 (inside the 
dashed lines), I//3///0I = 0.1 (solid lines), and 
l/^s/z^ol = 0-01 (dotted lines). 



A number of statistically less significant points appear outside the drawn regions in 
Fig. 5.18 and are not depicted. They correspond to points with m^-^ — nib — ^^.^ < 10 
GeV which appear above the diagonal line, and points with cos9^ < 0.1 which appear 
below the horizontal line corresponding to the lowest values of tan/5. In the last case, 
our approximation in eq. (5.17) does not work any more. On the other hand, eq. (5.17) 
predicts very well the behavior of r(r)/r(x^) if cos^^^ > 0.1. For example for f^s/fio = 1, 
or equivalently sin,^' = I/V2, we expect from eq. (5.17) a maximum value of order 1 for 
large tan/5 {hb ~ 1) and a minimum value of order 10~^ for small tan/5 {hf ^ 10"'^), 
and this is confirmed by Fig. 5.18. High values of the L violating branching ratio for 
large /i^ values are highly restricted for large tan/5. This can be understood as follows. 
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In the case of /xaZ/io = 1 and tan/? = 60 acceptable neutrino masses are obtained only 
if sin^ ~ 1. On the other hand, in this regime we find from eq. (1.26) that the AB 
term is large because of the high value of tan f3, and that the Am^ term is large because 
becomes negative and Am^ = mj^^^ — Ml grows in magnitude. This way, acceptable 
neutrino masses are achieved only with cancellation within more than one decade. In any 
case, we think that Fig. 5.18 is very conservative considering that in SSSM-SUGRA with 
unification of top-bottom-tau Yukawa couplings, the large value of tan/3 implies that a 
cancellation of four decades among vev's is needed. 

The width of the band in Fig. 5.18 reflects the degree of correlation between the 
ratio T(ti 6r)/r(ti bxf) and the neutrino mass under the mentioned conditions. 
Note that one would have an indirect measurement of the neutrino mass if this ratio 
were determined independently. The band will open to the left if one allows a stronger 
cancellation between the terms in AB and Am^. On the other hand it will open to the 
right if the universality between Bq and B^ is relaxed. This is shown in Fig. 5.19 where 
we plot the ratio r(r)/r(x^) in RII as a function of the z/3 mass, but without imposing 
universality between B^ and Bq. If we accept cancellation within one decade between the 
AB and Am^ terms, then the allowed region is at the right and below the corresponding 
dashed tilted line. If a larger degree of cancellation is accepted, the left boundary of 
the allowed region moves to the left as indicated in the figure, enhancing the L violating 
channel. In addition if we accept only a decade of cancellation between the two terms that 
contribute to the 1/3 mass, then our approximate formula which predicts the minimum 1/3 
mass in eq. (1.29) works very well. 




-1000-800 -600 -400 -200 200 400 600 800 1000 

Ao (GeV) 

Figure 5.20: Universality condition Bq = at the unification scale as a function of Aq. As 
tan P increases, the allowed values of Aq are more constrained. 

In summary, in this subsection we have shown that even in region II, where the L 
conserving decay mode ti bxt is also open, the L violating decay mode ti ^ br can 
be comparable to ii ^ bxt for large tan (3 and fi^, and relatively close to the chargino 
production threshold. In general, this implies a large neutrino mass unless a cancellation 
is accepted between the two terms contributing to the tree level neutrino mass. In ad- 
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dition, the non-universality of the Bq and terms at the GUT scale does not increase 
appreciably the allowed parameter space, except at large tan/?. The main consequence 
of this non-universality is to restrict the allowed values of at large tan /?. In the next 
subsection we study the effects introduced by the non-universality of and 

Effects of non— Universality 

We now study the effect of possible non-universality of soft-breaking SUSY parameters 
on our previous results. In particular, the non-universality between mj^^ and at the 
GUT scale. The SUGRA model with universality at the unification scale, while highly 
predictive, rests upon a number of simplifying assumptions which do not necessarily hold 
in specific models due to the possible evolution of the physical parameters in the range 
from Mpianck to Mqut- Specifically, there are several models in the literature with non- 
universal soft SUSY breaking mass parameters at high scales. A recent survey can be 
found in [193] , where several models such as based on string theory, M-theory, and anomaly 
mediated supersymmetry are analyzed. For this reason we find interesting to explore here 
the effects of non-universal soft terms. 

The SUGRA spectra are typically found for given values of mi/2, itiq, Aq, tan/? and 
Sgn(/io). In our case we have in addition sin^' (or equivalently, /is). The value of is 
determined by the previous parameters through the minimization conditions. In addition, 
a relation between Aq and the ratio B^/Bq at the GUT scale (which indicates the degree 
of universality) emerges. This relation can be seen in Fig. 5.20 for /^s/fJ^o = 1 and the 
values tan/? = 3, 40, and 60, for ml^^ = M|. The relation becomes more restrictive as 
tan/? is increased, starting from —1000 < Aq < 1000 GeV allowed for tan/? = 3, to 
a single value compatible with unification for tan/? — 60. Another way to enhance 
the L violating channel, enlarging the band toward the left in Fig. 5.18, is by relaxing 
the universality between m^^^ and M| at the GUT scale. In Fig. 5.21 we plot the ratio 
ml^yMl at the weak scale as a function of the same ratio at the unification scale Mqut for 
tan/? = 3. The shaded region is allowed, implying a maximum value for the ratio Tn\^/M'j^ 
at the weak scale for a given value of the ratio at the GUT scale. We see from Fig. 5.21 
that a relaxation of universality of 0.5% or more is enough to make {m\^/Ml)yjeak — 1 
possible, meaning that smaller neutrino masses are attainable without having to rely on 
a cancellation between the Am^ and /\B terms or small values of T{ti br). 

However as we increase tan/? the maximum value of m\^/Ml decreases, and thus, 
the required non- universality between rriLo and at unification scale grows drastically. 
In Fig. 5.22 we show the ratio Tn\^/Ml at the weak scale as a function of tan/?. We 
appreciate clearly the growing of |Am^|^j„ with tan/?. We remind the reader that this 
kind of non-universality in the soft terms is not uncommon in string models [194], or 
GUT models based on SU{5) [195] or 5*0(10) [196] for example. There are in fact some 
5'O(10) models for non- universality of the GUT scale scalar masses which naturally favor 
fight neutrino mass [196]. 

The effect of non-universality it is also explored in Fig. 5.23 where it is shown the 
relation between the neutrino mass and the parameter sin^ for Hs/fio = 1. Two different 
bands are shown: one for tan/? = 3 and T{t)/T{xi) = 2 x 10~^, and a second one for 
tan/? = 46 and r(T)/r(x^) = 0.4 ± 0.2. The required degree of universafity at the GUT 
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Figure 5.21: Comparison between the ratio 
rnj^^jM^ at the weak and the unification scales 
for tan f3 = 3. Universahty at the unification 
scale, mj^^/M'^ = 1, implies a maximum value 
for this ratio at the weak scale. 



Figure 5.22: {m\^jMl) evaluated at the weak 
scale as a function of tan (j. This ratio is al- 
ways less than one and decreases with tan (3. 




Figure 5.23: Minimum value of the 1^3 mass as a function of sin^ for different values ol rriLo/ML 
at the GUT scale and two values of tan /5. The ratio IJ-3/ fio is fixed to the indicated value, leading 
to a nearly constant value for T{ii 6r)/r(ti — > bxl). Here we assume that the two terms 
contributing to the fs mass cancel to within an order of magnitude. 
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scale is indicated inside the bands. For example, in order to have neutrino masses of 
the order of eV for tan/? = 3, m|y needs to be at least 0.2% larger than M|. Similarly, 
for tan/3 = 46 we need a m^^ twice as large as M| at the GUT scale in order to have 
neutrino masses of 1 eV. We stress the fact that for Fig. 5.23 we have conservatively 
accepted cancellation at the level of one order-of-magnitude only. 

In summary, the lesson to learn here is that non-universal soft SUSY breaking terms at 
the GUT scale have the potential of making it easier to reconcile sizeable L violating effects 
in the stop sector with very small neutrino masses, without resorting to cancellations. 



5.2 Three Body Decays of stop 

We discuss the phenomenology of the lightest stop in the SSSM. In this class of models 
we consider scenarios where the L violating two-body decay ii ^ h competes with the 
leading three-body decays such as ti — > W^h'j(^. We demonstrate that the L violating 
decay can be sizable and in some parts of the parameter space even the dominant one. 
Moreover we discuss the expectations for fi — > 6 and ii — > e"*" h. The recent results from 
solar and atmospheric neutrinos suggest that these are as important as the r^h mode. 
The ti ^ h decays are of particular interest for hadron colliders, as they may allow 
a full mass reconstruction of the lighter stop. Moreover these decay modes allow cross 
checks on the neutrino mixing angle involved in the solar neutrino puzzle complementary 
to those possible using neutralino decays [167]. For the so-called small mixing angle 
or SMA solution ti ^ h should be negligible, while for the large mixing angle type 
solutions all ti h decays should have comparable magnitude. We first explore the 

extent to which the decay ti — > 6t can be sizeable when compared with the 3-body 
decay modes. Moreover, we discuss the connections between the decay modes fi — > 6 Z+ 
and neutrino physics, in particular we discuss a possible test of the solution to the solar 
neutrino puzzle. The appendixes in [106] contain complete formulas for the total widths 
of the three-body decay modes as well as for the couplings. 

Notice that in SSUGRA models [31] mj,^ is calculable through the RGE evolution and 
one finds in this case cancellations up to two orders of magnitude for the combination 
A3 = —fisVo + /J'qVs. In general the smallness of rrii,^ requires relatively small 113 as in 
ref. [50]. The remaining two neutrinos acquire mass radiatively. Rigorous quantitative 
results were given in the second paper in ref. [34] . Typically they are hierarchically lighter 
than the heaviest neutrino, whose mass arises at the tree-level. This way one accounts 
for the observed hierarchy between the solar and the atmospheric neutrino mass scales. 

In this section we present our numerical results for the branching ratios of the lighter 
stop ^1. Here we consider scenarios where all two-body decays induced at tree-level are 
kinematically forbidden except the b decays. Before going into detail it is useful to have 
some approximate formulas at hand [31]: 

Ffti ^ 6t) fti !v ; — (5.18) 

^ ^ 167r 167r|/xop 



mlo 



^ cx",) ^ Fht{6^2Cosei- SAsmeiYflmi^ [l-^j , (5.19) 
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where F=^ {\og{mGVT/mz)K,hKt,/167r^f ~ 6 x 10"^ h = V2{teinew Nn + m2)/G 
and the parameter (5^2 is given by 

For the minimal SUGRA models one finds d„^2 = 0{1) which is basically independent 
of the initial conditions due to the ttiq dependence both in the numerator and in the 
denominator. Finally we have 

mt{Ah + I At) 

"A — 2 2 

m~ — mf 

The complete formulas are given in [31, 70], while for the three-body decays they are 
given in Section 4.2.2 . They reduce to the ones given in [73, 74] for vanishing L violating 
parameters. 

We have fixed the parameters as in [73] to avoid color breaking minima, while in 
the top squark sector we have used m^^, cos^f, tan (3, and fiQ as input parameters. For 
the sbottom sector we have fixed Mq^Mj^ and as input parameters whereas for the 
charged scalars we took mpo, M^,Mi, and Aj. as input ^. In addition we have chosen 
the L violating parameters /i^ and in such a way that the heaviest neutrino mass is 
fixed with the help of eq. (1-24). For simphcity, we have also assumed that the soft SUSY 
breaking parameters are equal for all generations. 

In order to get a feeling for minimal size of branching ratios that can be measured let 
us first shortly discuss the expected size for the direct production of light stops at future 
coUiders. One expects for example at the LHC a production cross section of ~ 35 pb 
for 220 GeV stop mass. Therefore, once the full luminosity has been reached, one has to 
expect approximately 3.5 10^ events per year. The corresponding stop production cross 
section at a future e+e^ linear collider of 800 c.m.s. energy is of (9(10 — lOOfb) [78]. For 
an integrated luminosity of 500 fb~^ per year one can expect 0(10"^) stop pairs per year. 
This implies that branching ratio as low as 10~^ can in principle be measured. 

We consider first the simplest case of one generation model which, as already men- 
tioned, is sufficient to describe the relative importance of the — > r"*" 6 decay mode 
relative to the possible 3-body decay modes 

ii ^ W+bxl 

ii - SthxH 

h &'5^T+, (5.20) 

ii — > bi>i l'^ {I — e, ii) . 



^nipo plays here the same role as m^o in the MSSM case [73]. 



5 Numerical Results 



73 



In general the important final states are those that conserve L. For example, for the 
case of decays involving the most important are those in which the scalars will be 
mainly a stau. Due to the fact that the existing bounds on the MSSM sneutrinos are 
below 100 GeV there exists the possibihty that the sneutrino has nearly the same mass 
as one of the Higgs boson. Similarly it could be that the charged MSSM boson has nearly 
the same mass as one of the staus. This implies large mixing effects even for small L 
violating parameters [84, 85, 105]. We therefore have used the complete formulas for the 
3-body decay modes which are presented in the Appendix of [106]. The latter include L 
violating decays such as ii — > W'^bu^. In addition to the above mentioned decays there 
is also ii — > bZ^T'^. This decay mode is kinematically suppressed compared to ii ^ br'^ 
and there is no possible enhancement due to a mixing with an L conserving final state. 
Therefore it can be safely neglected. 

In Fig. 5.24 we show the branching ratios for the ti as a function of cos% in different 
scenarios. In order to calculate the partial width for the decay ii — > cXi we have taken 
the formula given in ref. [71]. According to the analysis performed in [31], where the full 
calculation was done in the SUGRA scenario with universality at the unification scale, 
the result obtained with the present approximation should be taken as one upper bound. 
This implies also that the shown branching ratio for ti — > 6r"^ can be viewed as a lower 
bound. The parameters and physical quantities used in Fig. 5.24 are given in Tab. 5.3. 
For the case of Fig. 5.24(a) we have fixed in addition the L violating parameters such that 
rrii, ^ 1 eV. With this choice of parameters S^, S^, P3, and are mainly the MSSM 
Higgs-bosons whereas S^, P2, S2 , and are mainly the MSSM sleptons of the third 
generation. In the plot we show the various branching ratios of the lighter stop summing 
up those branching ratios for the decays into sleptons that give the same final state, for 
example: 

ii^bu^el 6 e+ z^e Xi , h^be^ ^ b e+ i^e Xi ■ 

The branching ratios for decays into fli or i/^ are practically the same as those into 
e-L or i>g. Note that the energy spectrum of the leptons in the final will be somewhat 
different depending on whether the scalar in the intermediate step is charged or neutral. 
This offers in principle the possibility of determining the branching ratios of the different 
decay chains even if the final state topology is common. Note, that states containing 
scalars or neutralinos will lead to additional jet and/or lepton multiplicities absent in the 
MSSM. 

In Fig. 5.24(b) the slepton mass parameters are chosen such that decays into scalars 
are kinematically forbidden. Here we display the channels ti — > bW'^Xi: h ~^ ^'^^ 
and tl — > cx5- The remaining modes, such as ti — > fcS'J'r"'", turn out to be completely 
negligible. In both cases, with and without sleptons in the final state, one can see that 
in general the three body mode ii — > bW^'j(^, dominates except for a somewhat narrow 
range of negative cos%. However, the branching ratio for ii — > 6r+ is above 0.1% for 
most values of | cos 6i\ implying the observability of this mode. Most importantly, note 
that even in the parameter ranges where the three-body decay mode is dominant, its 
resulting signature is rather different from that of the MSSM due to the fact the lightest 
neutralino decays into SM-fermions, leading to enhanced jet and/or lepton multiplicities, 
as discussed in detail in [103, 167]. In the remaining part of this section we assume that 
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Figure 5.24: Branching ratios for the ti as a function of cos Of for different scenarios. 
We have fixed in a) m^^ = 1 eV, b) m^^ = 1 eV, > 225 GeV, > 225 GeV, c) 
m^3 = 0.06 eV, > 225 GeV, > 225 GeV, d) Branching ratios for the ti as a 
function of cos^^- for tan/3 = 3. m^.^ = 0.06 eV, > 225 GeV, > 225 GeV. All the 
other inputs are given in Table 5.3. 
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Table 5.3: Input parameters and resulting quantities used in Fig. 5.24. 

3-body decays into scalars are kinematically forbidden. 

In Fig. 5.24(c) the L violating parameters are fixed in such a way that the heaviest 
neutrino mass is in the range suggested by the oscillation interpretation of the atmospheric 
neutrino anomaly [26]. 

B( fj ^ b T ) 




5 10 15 20 25 30 



Figure 5.25: Branching ratios for ti decays for m^-^ = 220 GeV, |Uo = 500 GeV, M2 = 240 GeV, 
and rriy = 100, 1 and 0.06 eV. The branching ratios are shown as a function of tan/3, (cos^^ = 
-0.8) 

In Fig. 5.24(d) we show the same scenario as in Fig. 5.24(c) but for tan/3 = 3. The 
branching ratio into br now increases, whereas the branching ratio into cXi decreases. 
This is easily understood by inspecting eqs. (5.18) and (5.19). Indeed for the br case 
the partial width is proportional to hi, whereas for cXi it is proportional to h^. This 
implies that the partial width for ii — > cXi grows faster with tan P than the width for 
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ti ^ br. This is also demonstrated in Fig. 5.25 where we show the tan/? dependence of 
the branching ratio for the decay of ti into 6r+ for several values of the neutrino mass. 
For m^^ = 0.06 eV the BR{ii — > 6t) is still above 0.1% if tan /3 is not too large, as favored 
by the explanation of the neutrino anomalies in this model [34] . As seen from the figure, 
the the ti — > fer"*" branching ratio is also somewhat correlated to the mass. Should one 
add a sterile neutrino to the model [59], then the neutrino state could in principle be 
heavier than assumed above, favoring ti ^ r+ 6 decay mode. 

Let us now turn to the general three neutrinos case. There are new features that arise 
in this case, as opposed to the 1-generation case considered so far. In this model the 
solution to the present neutrino anomalies implies that all the /li are of the same order of 
magnitude [34]. 

Two further important results of [34] are that the atmospheric neutrino angle is con- 
trolled by the ratio {—112^0 + l^oV2) / {— I^^Vq + HqV^) and that the solar mixing angle is 
controlled by (//i//^2)^- One can get approximate formulas for the decay widths fi — > fee"*" 
and ii ^ h 11^ similar to eq. (5.18) by replacing /X3 by /Xi^2- This implies that (i) The 
decays into fee"*" and 6//+ are as important as the decay into 6r+. (ii) The decays 
ti 6e+ and ti b fi^ are related with the solar mixing angle. Moreover, we find 
that Xli=e ^ T -'^(^1 ~^ bl^) in the 3-generation model is nearly equal to r(ti br^) in the 
1-generation model provided that Y^i=i lA identified to /Xg in the 1-generation model. 

In Fig. 5.26 we show the ratio of BR(ii — > be'^)/BR{ii — > 6/^"*") versus {jii/ IJL2Y 
for different values of cos 9];. For definiteness we have fixed the heaviest neutrino mass 
at the best-fit value indicated by the atmospheric neutrino anomaly. One can see that 
the dependence is nearly linear even for rather small cos Of:. For |cos^^j| < 10~^ the 
approximation in eq. (5.18) breaks down and additional pieces dependent on sin^^^ [31, 
70] become important, leading to the non-linear dependence. One sees from the figure 
that, as long as cos^^^ > 10^^ there is a good degree of correlation between the branching 
ratios into BR(ti 6e+) and BR(ti bfi^) and the ratio (yUi//X2)^- Thus by measuring 
these branchings one will get information on the solar neutrino mixing, since tan^ 6soi 
is proportional to {^i/ 1^2)^ [34] which makes it a rather important quantity. For the 
so-called small mixing angle or SMA solution of the solar neutrino problem we expect 
ti ^ b to be negligible. In contrast, for the large mixing angle type solutions (LMA, 
LOW and QVAC, see ref. [26] and references therein) we expect all ti b decays to 

have comparable rates. As a result in this model one can directly test the solution to the 
solar neutrino problem against the lighter stop decay pattern. This is also complementary 
to the case of neutralino decays considered in [167]. In that case the sensitivity is mainly 
to atmospheric mixing, as opposed to solar mixing. Testing the latter in neutralino decays 
at a collider experiment requires more detailed information on the complete spectrum to 
test the solar angle [167]. In contrast we have obtained here a rather neat connection of 
stop decays with the solar neutrino physics. 

Note, that this result is much more general than the scenarios discussed in this work. 
It is of particular importance in scenarios where only the L violating decays and the 
decay into are present [31, 70]. Similarly, the other ratios of the final states bl~^ are 
proportional to the square of the ratio of corresponding /Xj provided that cos 9i is not too 
small. 
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Figure 5.26: Ratio of branching ratios: BR{ti be~^)/BR{ti b^~^) as a function of 

{1^1/^2^ for TTii^ = 220 GeV, /xq = 500 GeV, M2 = 240 GeV; Icos^^l > 0.1,0.01,10-3, 
m^^ = 0.6 eV. 
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Chapter 6 
Conclusions 



We have seen in this work that the SSSM is a well motived model which also provides a 
sohd explanation of the neutrino anomahes [34, 38] which can be tested at accelerators 
experiments. It reproduces the same signals that general Sp-SSM but with a minimum set 
of L violating terms and automatically satisfies all the experimental bounds on induced 
dimension-4 L violating operators. We have explicitly studied the phenomenology of 
neutralino and stop decays. 

Concerning the Neutralinos, we have studied the production of the lightest neutralino 
x5 at LEP2 and the resulting phenomenology in models where an effective superrenor- 
malizable term parameterizes the explicit breaking of L. We have considered supergravity 
scenarios which can be explored at LEP2 in which the lightest neutralino is also the light- 
est supersymmetric particle. We have presented a detailed study of the LSP x5 decay 
properties and studied the general features of the corresponding signals expected at LEP2. 
A detailed investigation of the possible dctcctability of the signals discussed in Tab. 5.1 
taking into account realistic detector features is beyond the scope of this thesis. Clearly, 
existing LEP2 data are already probing the part of the parameter region which corre- 
sponds to approximately m^o < 40 GeV. Finally, we note that, in addition to important 
modifications in the x? decay properties, the SSSM lead also to new interesting features in 
other decays, such as charged [85] and neutral [84] Higgs boson and slepton decays, stop 
decays [31, 70], and gluino cascade decays [77, 197]. In addition wc have shown that the 
SSSM can not be confused with gauge mediated supersymmctry breaking and conserved 
R-parity due to the absence of several final states in the GMSB case. 

Concerning the stops we have studied the two-body decays of the lightest top squark in 
SUGRA models with and without L violating terms. We have improved the calculation for 
the decay ti — c x° by numerically solving the renormalization group equations (RGE's) 
of the SSSM including full generation mixing in the RGE's for Yukawa couplings as well 
as soft SUSY breaking parameters. The decay-width is in general one order of magnitude 
smaller than the one obtained in the usual one-step approximation. This result will 
therefore enlarge the regions of parameter space where the four-body decays of lightest 
stop dominate over the decay into a charm quark and the lightest neutralino. As a result 
it will affect the present experimental lower bound on the ti mass even in the L conserving 
case [76] . In the SSSM of course new decay modes appear and, as we have shown, they 
can be sizeable. In fact we have shown that the lightest stop can be the LSP, decaying 
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with 100% rate into a bottom quark and a tau lepton. We have shown that the decay 
mode ti ^ br dominates over ti — > c even for neutrino masses in the range suggested 
by the simplest oscillation interpretation of the Super-Kamiokande atmospheric neutrino 
data. This result would have a strong impact on the top squark search strategies at 
LEP [198] and TEVATRON [199], where it is usually assumed that the ^ cx° decay 
mode is the main channel. In addition to the signal of two jets and two taus present when 
the two produced stops decays through the L violating channel, one expects a plethora 
of exotic high-multiplicity fermion events arising from neutralino decay, since such decay 
can happen inside the detector even for the small neutrino masses in the range suggested 
by the z/^ to Ur oscillation interpretation of the atmospheric neutrino anomaly [34] . 

In addition the L violating decay ti ^ hr^ may compete with three-body decays. 
Specifically, we have found that for m^^ < 250 GeV there are regions of parameter where 
ti — > 6 r+ is an important decay mode if not the most important one. This implies that 
there exists the possibility of full stop mass reconstruction from t'^ t~ hh final states, 
favoring the prospects for its discovery. In contrast, in the MSSM the discovery of the 
lightest stop might not be possible at the LHC within this mass range. This implies 
that it is important to take into account this new decay mode when designing the stop 
search strategies at a future e'^e" Linear Collider. The SSSM also imply additional leptons 
and/or jets in stop cascade decays. Looking at the three generation model the decays into 
^1 — > imply the possibility of probing and thus the solar mixing angle. This 

complements information which can be obtained using neutralino decays. In the latter case 
the sensitivity is mainly to the atmospheric mixing, as opposed to solar mixing [167]. In 
SSSM neutralino decays is ideal to test the atmospheric anomaly at a collider experiment, 
while stop decays provide neat complementary information on the solar mixing angle. 
Obtaining solar mixing information from neutralino decays would require more detailed 
knowledge on the supersymmetric spectrum, since it would be involved in the relevant 
loop calculations of the solar neutrino mass scale and mixing angle. By combining the 
two one can probe the parameters associated with both solar and atmospheric neutrino 
anomalies at collider experiments. 



Appendix A 

Constraints on horizontal charges 



In this Appendix we derive the general expressions for the individual field charges sat- 
isfying the set of 13 phcnomenological and theoretical constraints corresponding to the 
six mass ratios for the quarks and the charged leptons plus the two quark mixing angles 
summarized in Table 2.3; the two relations provided by the absolute value of the masses 
of the third generation fermions; one phenomenological assumption about the charge of 
the fiQ term (2.20); one theoretical constraint corresponding to the consistency conditions 
for the coefficients of the mixed linear anomalies and one additional constraint from the 
vanishing of the mixed anomaly quadratic in the horizontal charges. As discussed in sec- 
tion 2.3.1, this leaves us with four free parameters that we choose to be {i = 1,2,3) 
and X. We obtain 



15 (7 + x) 



180 - 45a; - + Qi3(41 + 5a;) - 7L23 + L: 



23 



+ni(2 + X + L23) + ^2(9 + x- L23) + ^3(9 + x] 



15(7 + 2;) 



20 + 50a; + 6a;^ + 18^13 - 2IL23 + 3L 



23 



-ni(29 + 2a; - 3L23) - 712(8 + 2a; + 3L23) + n3(97 + 13a;) 



where L23 = L23 -|- ^23 and Q13 parametrized the two different possibilities for the quark 
and lepton charge differences given in Tab. (2.3). In terms of Qs and L3 and of our four 
free parameters we have 



Hu = ris- L3 
Lq = — 1 -I- Hu 



U3 = -Q3 - 

ds = -Q3 + x- Lq 

h = -L3 + X - Lq 



and from these all the other individual charges can be straightforwardly determined from 
the charge differences in Tab. (2.3). The solution for the charges in model MQl+MLl for 
the preferred values rii = n2 = ns = —8 and a; = 1 is given in Table A.l. For the no = 
case and the same value for the neutrino mass which implies in this case rii — —6 — x, we 
present in Table A. 2 the individual charges in in model MQl-l-MLl for a; = 
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Qi 


Q2 


Qs 


Ul 


U2 


"3 


di 


d2 




Ll 


L2 




h 


h 


^3 


Hu 


Lo 


161 


131 


71 


103 


58 


28 


18 


23 


23 


113 


113 


113 


98 


53 


23 


127 


97 


30 


30 


30 


15 


15 


15 


5 


5 


5 


30 


30 


30 


15 


15 


15 


30 


30 



Table A.l: The anomaly free set of charges of model MQl+MLl for x — 1 and rii — n2 — 
n3 = -8 
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23 


18 
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28 


13 
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14 
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19 


19 


28 


18 
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11 
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5 


5 
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5 


5 


5 


5 


5 


5 


5 


5 


5 


5 


5 


5 


5 



Table A. 2: The anomaly free set of charges of model MQl+MLl for uq — 0, x — and 
rii — n2 — — —6 



Appendix B 

Complete formulas for 2 body decays 



In this Appendix we derive the Feynman rules Ffqiq^ (involving a neutralino/tau- neutrino, 



a quark, and a squark) and Fj qiq'/^ (involving a chargino/tau, a quark, and a squark of 
different electric charge) in the case of three generations and L violating in the third 
generation. This is a generalization of the Feynman rules contained in [200] , which hold 
for the MSSM and for one generation of quarks and squarks. 



B.l An explicit calculation 

Before the general case, we will calculate here ti — > x? 

From eq. (1.4), the lagrangian in terms of the two body decays are 



1 ~, 



-=tl [9{iX^)tL + 9'ht{iX')tL + h.c] + -^g'ht [{iXyJl + tLi-iX')^] 



V2 



V2' 



In terms of the four component spinors 



-ht{i*j,tLH'^ + tLtRH^)} (B.l) 



B 



iX' 



-iX^ 
iX' 



H„ 



~0 



(B.2) 



we have 



gWM + g'htBPLt 



httPLBtR + KBPRti*ji 

grrit 



\/2mw sin (5 



HuPLtf^ + iPLH^tL (B.3) 



The result in the mass basis is finally 

= 9 [t [hljPL + fljPR] x'ii + h.c} 



(B.4) 
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where 



+ 



sin 6 1 



-2V2 . 
sm 



+ 



{sm(3N'j^-cos(3N*^) 



cos 9 



w 



cos9^ 



sin 9^ 



where we have redefined the neutrahno diagonahzation matrix as 



cos 9w sin 6'vi/ 
— sin 9w cos 6*14/ 



J2. 







COS — sin /5 










sin /9 cos /5 




_iV,4_ 



We now calculate the two body width decay T(ti ^ t + xf) 



r(ti^t + xi) 



ti 



where A(a, b,c) = {a — b — c)^ — 4ac. The corresponding Feynman diagram is 




Consequently we can write down the amplitude as 

M=g [h\^u{pi)PLv{p2) + fl,u{pi)PRv{p2)] 

The final result is 



B.2 Couplings Neutralino— Squark-Quark 



(B.5) 



(B.6) 



[B.7) 



(B.8) 



(B.9) 



(B.IO) 



Following [201] we work in a quark interaction basis where dL,R = r, ul = Ku^^^, and 
ur = M^j. (we denote q and g° the mass and current eigenstates respectively), as opposed 
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to Ref. [202] where a more general basis is used. In addition, we implement the notation 

ql pi = ^ for the interaction basis. 

The starting point is the following piece of the Lagrangian 



■'uuFO — —gu 



sin Ow^uN'ji + 



cos 6w ^ ^ 



i -et/sin^ ew)N'j2 



UiL 



m. 



ou 

ij 



\/2mw sin (3 sin 9 



sinOweuN'j^ + (-et/sin^ Ow)N'*j2 



om 



\^mw sin P sin 9 



cos 9w 
N'*j,ujAPLF'j + h.c. 



UiR 



(B.ll) 



written in the quark interaction basis. The 5x5 matrix N' diagonalizes the neu- 
tralino/ncutrino mass matrix in the (7, Z, Lg, z/^) basis as defined in [85], with the 
index J = 1...5. The 3x3 up-type quark mass matrix m"^ is not diagonal, with the 
indexes i,j — 1, 2, 3. 

In order to write the above Lagrangian with mass eigenstates we use the basic relations 
mentioned before, in particular, u^j^ — {K^Y^UjL, which implies that ufj^ — Uj^K^^ We 
need the following relations: 



ki 



Uk 



(B.12) 



*^U\ki~ 



where i,j — 1,2,3 label the quark fiavors, k = 1...6 labels the squarks, and mF = 
diag {rriu, rric, rrit} is the diagonal up-type quark mass matrix. In this way, the Lagrangian 
in eq. (B.ll) can be written as 



where the different couplings are 



(B.13) 



fyjki 
'^OUL 



f-tjki 
^OUR 



sin 9w&uN'*i + 



cos 9—^'^ 



a-eusin^9w)N'* 



i2 



ki 



sin 9weuN'ji + 



cos 9 



w 



(-ec/sin^ 9w)N'j2 



■^ki 
UR 



(B.14) 



^ouR = N'*Jj:uRhuT'' 



and hu =diag (m„, mc, mt)/(\/2mi^ sin /5 sin 6*). Graphically, the FjUiUk Feynman rules 
are given by 
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The analogous Feynman rules in the MSSM are obtained by replacing by 
interpreting the matrix as the usual 4x4 neutralino mixing matrix, and by setting 
= 7r/2 in the formula for the Yukawa couplings. 

Similarly, replacing all u{u) by d{d) in eq. (B.ll) and starting from 



^qq'F+ —gdi 



out 



^/2mw sin /? sin 9 



+ gdi 



^/2mw cos l3 sin 9 

D] 

^/2mw cos P sin 9 



Vj2UjR - VjiUiL 



PrPj 



Uj2djR — UjidiL 



PrF} 



my 



\/2mw cos P sin 9 



VhhLPLFj + h.c 



(B.15) 



we can obtain the complete Feynman rules for the neutralino/tau-neutrino and chargino/tau 
with quarks and squarks. The results, that complements the obtained in [201], are 

Neutralino-(d)quark-(d)squark 




M{V2G*^dI + H*^dr)Pr - iV2G*^DR - H*^dl)Pl] 
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The mixing matrices God and Hqd are defined as 



sin 9weDN'*^ H b~(^3D - sin^ 9w)N'*2 

■' cos Ow 



DL 



cos 6^ 



-CD sin^ Ow)N' j2 



r 



ki 
DR 



(B.16) 



B.3 Couplings Chargino-Squark-Quark 

Here we give the couphngs that were used in Sections B.2 and 4.2.2 
Chargino/tau-(d)quark-(u)squark 




-^gi-C-')[iG^^l 



HIjr)Pl - HiflPR] 




*jki - 



where C is the charge conjugation matrix (in spinor space) and the mixing matrices Gu 
and Hu are defined as 



ki 



ki 



(B.17) 



Chargino/tau-(u)quark-(d)squark 




jki ■■ 



jjki 



R\ 



^9[{Gt,'-HZ')PR-H*^rPL]C 



f-*jki 



*jki j 
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(B.18) 



where the mixing matrices Go and Hd are defined as 

In order to derive the decays widths we write, for example eq. (B.13) as 

jCuufo = gu^rJ^'PR + ht'PL)F^Uk + h.c (B.19) 
The result is presented in Sec. 4.2.1. 
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